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A Model of Turbulence in Magnetized Plasmas: Implications for 
the Dissipation Range in the Solar Wind 

G. G. Howes/ S. C. Cowley,^' ^ W. Borland, G. W. Hammett,^ E. Quataert/ and 
A. A. Schekochihin^ 

Abstract. This paper studies the turbulent cascade of magnetic energy in weakly col- 
lisional magnetized plasmas. A cascade model is presented, based on the assumptions 
of local nonlinear energy transfer in wavenumber space, critical balance between linear 
propagation and nonlinear interaction times, and the applicability of linear dissipation 
rates for the nonlincarly turbulent plasma. The model follows the nonlinear cascade of 
energy from the driving scale in the MHD regime, through the transition at the ion Lar- 
mor radius into the kinetic Alfven wave regime, in which the turbulence is dissipated 
by kinetic processes. The turbulent fluctuations remain at frequencies below the ion cy- 
clotron frequency due to the strong anisotropy of the turbulent fluctuations, fc|| k± 
(implied by critical balance). In this limit, the turbulence is optimally described by gy- 
rokinetics; it is shown that the gyrokinetic approximation is well satisfied for typical slow 
solar wind parameters. Wave phase velocity measurements are consistent with a kinetic 
Alfven wave cascade and not the onset of ion cyclotron damping. The conditions under 
which the gyrokinetic cascade reaches the ion cyclotron frequency are established. Cas- 
cade model solutions imply that coUisionless damping provides a natural explanation for 
the observed range of spectral indices in the dissipation range of the solar wind. The dis- 
sipation range spectrum is predicted to be an exponential fall off; the power-law behav- 
ior apparent in observations may be an artifact of limited instrumental sensitivity. The 
cascade model is motivated by a programme of gyrokinetic simulations of turbulence and 
particle heating in the solar wind. 



1. Introduction 

In the study of turbulence in magnetized plasmas, the so- 
lar wind provides a unique opportunity to characterize the 
nature of turbulent plasmas by making detailed dynamical 
measurements through in situ satellite observations. But us- 
ing the solar wind as a turbulence laboratory is complicated 
by the fact that the solar wind is only weakly coUisional, 
so a kinetic description is necessary [Marsch,, 1991; Marsch, 
2006]. A useful first step in connecting theoretical ideas 
to observational data over the range from the large energy- 
containing scales of the turbulence to the small scales where 
the turbulence is dissipated is to have a simple model com- 
bining our understanding of magnetized turbulence at both 
fluid and kinetic scales. This paper presents such a model. 

This paper follows Howes et al. [2006] and Schekochihin 
et al. [2007] as part of a program to study kinetic turbulence 
in astrophysical plasmas using nonlinear gyrokinetic simu- 
lations. Here we attempt to model the observed magnetic 
energy spectrum in the solar wind using a minimal num- 
ber of ingredients — namely, finite Larmor radius efi^ects and 
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kinetic damping via the Landau resonance. A lengthy eval- 
uation of the applicability of the gyrokinetic approximation 
to the solar wind plasma is presented along with a discussion 
of the relevant observational data. The gyrokinetic cascade 
model yields sensible results when compared to the obser- 
vational data, but much additional physics exists and may 
prove to be important; nonlinear simulations will play a crit- 
ical role in unraveling the complex interactions occurring in 
the turbulent solar wind. 

1.1. Turbulence in Magnetized Plasmas 

The ubiquity of turbulence in space and astrophysi- 
cal plasmas — e.g., the solar corona and solar wind, accre- 
tion disks around black holes, the interstellar medium — has 
driven a vast efi^ort to understand and characterize turbu- 
lence in a magnetized plasma. Progress in the understand- 
ing of magnetohydrodynamic (MHD) turbulence has accel- 
erated steadily over the decades. The pioneering work of 
Irostinilcov [1963] and Kraicfinan [1965] identified two key 
properties of MHD turbulence: first, even in a magnetized 
plasma with no mean field, the magnetic field of large- 
scale fiuctuations can behave efi'ectively as a mean field for 
smaller-scale fluctuations; and second, only oppositely di- 
rected Alfven wave packets interact nonlinearly. They as- 
sumed weak nonlinear interactions and an isotropic (with 
respect to the direction of the mean magnetic field) cascade 
of energy and derived an energy spectrum that scaled as 
Measurements of magnetized turbulence in labora- 
tory plasmas [Robinson and Rusbridge, 1971; Zweben et al., 
1979; Montgomery and Turner, 1981] and in the solar wind 
[Belcher and Davis, 1971] as well as the results of numeri- 
cal simulations [Shebalin et al., 1983; Oughton et al., 1994] 
brought attention to the anisotropy inherent in MHD tur- 
bulence in the presence of a mean field and inspired the 
early anisotropic theories of MHD turbulence [Montgomery 
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and Turner, 1981; Montgomery , 1982; Shebalin et al, 1983; 
Higdon, 1984]. Based on these ideas, Sridhar and Goldre- 
ich [1994] constructed a perturbative theory of weak MHD 
turbulence in which energy cascades only to higher perpen- 
dicular wavenumbers; but their neglect of the fcy — modes 
in resonant three- wave couplings provoked criticisms [Mont- 
gomery and Matthaeus, 1995; Ng and Bhattacharjee, 1996], 
and in response refined theories of weak MHD turbulence 
emerged [Goldreich and Sridhar, 1997; Ng and Bhattachar- 
jee, 1997; Galtier et al., 2000; Lithwick and Goldreich, 2003]. 

The anisotropic nature of weak turbulence implies that 
it generates fluctuations that become increasingly strong as 
the cascade proceeds to higher perpendicular wavenumber; 
combining this concept with the ideas of Higdon [1984], Gol- 
dreich and Sridhar [1995] proposed that a state of strong 
turbulence is eventually reached that maintains a critical 
balance between the (parallel) linear propagation and (per- 
pendicular) nonlinear interaction timescales. They pro- 
ceeded to develop a scaling theory of strong turbulence in 
incompressible MHD, hereafter referred to as the GS the- 
ory. Later studies have examined the effects of compressibil- 
ity [Lithwick and Goldreich, 2001; Cho and Lazarian, 2003; 
Vestuto et al., 2003], imbalanced Elsasser fluxes [Maron and 
Goldreich, 2001; Cho et al, 2002; Lithwick et al, 2007], and 
a host of other refinements [Oughton et al., 2004; Oughton 
and Matthaeus , 2005; Zhou and Matthaeus , 2005; Boldyrev, 
2005; Shaikh and Zank, 2006]. 

The GS theory for strong, incompressible MHD turbu- 
lence rests on two central assumptions: the locality of inter- 
actions in the wavenumber space [Kolmogorov, 1941], and 
the GS conjecture of critical balance. These imply a one- 
dimensional kinetic energy spectrum Ek{k±) oc k^^^^ (see 
§ 2.1). Using the consequent scaling for the nonlinear fre- 

2/3 

quency ujni ~ k±v± oc fcy and equating it with the lin- 
ear frequency for shear Alfven waves, uj — ±fc||t;A (critical 
balance), one finds fcy oc k'^^ . Therefore, the fluctuations 
become more anisotropic at high wavenumbers and the tur- 
bulence deep in the inertial range consists of nearly per- 
pendicular modes with k± ^ fcy. Numerical simulations of 
MHD turbulence with a dynamically strong mean field ap- 
pear to support these conclusions [Gho and Vishniac, 2000; 
Maron and Goldreich, 2001]. 

In real astrophysical plasmas, although the large scales at 
which the turbulence is driven may be adequately described 
by MHD, the turbulent fluctuations at the small-scale end 
of the inertial range often have parallel wavelengths smaller 
than the ion mean free path, fcyAmfpi ^ 1; therefore, a ki- 
netic description of the plasma is required to capture the 
turbulent dynamics (a systematic way to do this is given 
by Schekochihin et al. [2007]). It is, however, possible to 
show rigorously that Alfvenic fluctuations are essentially 
fluid in nature, satisfy the Reduced MHD (RMHD) equa- 
tions [Strauss, 1976], and remain undamped until their cas- 
cade reaches the perpendicular scale of the ion Larmor ra- 
dius, k±pi ~ 1 [Quataert, 1998; Lithwick and Goldreich, 
2001; Schekochihin et al., 2007]. In contrast, the compress- 
ible MHD fluctuations such as the slow, fast, and entropy 
modes require a kinetic description even at large scales and 
are damped both coUisionally by parallel ion viscosity [Bra- 
ginskii, 1965] at scales fcy Amfpi ~ 1 and coUisionlessly by the 
ion Landau damping [Barnes, 1966] at scales fcyAnifpi 3> 1. 
The cascade of these modes is not fully understood and they 
will not be considered here (see Schekochihin et al. [2007] for 
further discussion). 

Some fraction of the energy in Alfvenic fluctuations at 
k±pi ~ 1 is expected to launch a kinetic Alfven wave cas- 
cade to smaller scales [Gruzinov, 1998; Quataert and Gruzi- 
nov, 1999]. In the limit of anisotropic fluctuations with 
/C|| <C k±, this cascade, in the wavenumber range k±pi S> 1 



and fcxpe ^ 1, can again be described in fluid-like terms: 
the governing system of equations is known as Electron Re- 
duced MHD (ERMHD) [Schekochihin et al., 2007]. In the 
high-beta limit, it coincides with the anisotropic limit of the 
better known but nonrigorous fluid model called Electron 
MHD (EMHD) [Kmgsep et al., 1990], where the ions are 
assumed motionless and the current is carried entirely by 
the electrons. Assuming locality of interactions and a crit- 
ical balance between the linear and nonlinear timescales is 
again a plausible route towards a scaling theory of the ki- 
netic Alfven wave cascade. The outcome of this argument 
(explained in more detail in § 2.1) is the one- dimensional 
magnetic energy spectrum i5s(fcx) oc k^^'^ and a predicted 
wavenumber anisotropy that scales as fcy oc fc^'''^ [Biskmap 
et al., 1999; Gho and Lazarian, 2004; Krishan and Maha- 
jan, 2004; Shaikh and Zank, 2005; Schekochihin et al., 2007]. 
Simulations of EMHD turbulence support these predictions 
[Biskmap et al, 1999; Cho and Lazarian, 2004; Shaikh and 
Zank, 2005]. 

The MHD Alfven wave and kinetic Alfven wave cascades 
correspond to a nonlinear flow of energy to higher wavenum- 
bers across the {k±,k^^) plane, as illustrated in Figure 1. 
Imposing critical balance constrains the flow of energy to a 
one-dimensional path on the plane defined by a; = ujni, de- 
lineated on the figure as a solid line. Critical balance effec- 
tively defines a relation between k±, fcy, and the fluctuation 
amplitude 5B±; for strong turbulence [Goldreich and Srid- 
har, 1995], choosing values for two of these quantities spec- 
ifies the third. The freedom in this relation is conveniently 
parameterized by the wavenumber at which the fluctuations 
are isotropic, k± — fcy = ko. In Figure 1, adjusting the value 
of the isotropic wavenumber fcopi, or equivalently adjusting 
the amplitude of the fluctuations at a given value of k±pi, 
shifts the curve defining critical balance vertically; in this 
figure, we have chosen kopi — 10^*, which is the estimated 
value for the slow solar wind (see § 3.1). Since most treat- 
ments of magnetized turbulence assume an isotropic driving 
mechanism, we call the parameter kopi the isotropic driving 
wavenumber.^ We note here that critical balance, lj = ujni, is 
not strictly obeyed in MHD turbulence. Numerical simula- 
tions [Cho and Vishniac, 2000; Maron and Goldreich, 2001; 
Cho et al, 2002; Cho and Lazarian, 2003; Oughton et al., 
2004] suggest that the turbulent fluctuation energy does not 
simply reside along the path specified by critical balance 
but approximately fills the shaded region in Figure 1, where 
1^ ^ i-^ni- The one-dimensional cascade model presented in 
this paper can be thought of as integrating vertically over 
all at each value of k±, and then treating that energy as 
if it were concentrated on the path of critical balance. 

Although the fiuid MHD and EMHD treatments provide 
important insights into turbulence in magnetized plasmas, 
modeling accurately the nonlinear dynamics, in particular 
the coUisionless damping of turbulent fluctuations, requires 
a kinetic treatment. The inherently anisotropic nature of 
magnetized turbulence lends itself to description by an ap- 
proximation called gyrokinetics [Rutherford and Frieman, 
1968; Taylor and Hastie, 1968; Catto, 1978; Antonsen and 
Lane, 1980; Catto et al., 1981; Frieman and Chen, 1982; 
Dubm et al, 1983; Hahm et al., 1988; Brizard, 1992; Howes 
et al, 2006; Brizard and Hahm, 2007; Schekochihin et al., 
2007]. Gyrokinetics is a low- frequency expansion of kinetic 
theory valid in the limit of frequencies small compared to the 
ion cyclotron frequency ui <^ Qi, ion Larmor radius small 
compared to the length scale of the equilibrium quantities 
pi <^ L, wavevectors of electromagnetic fluctuations nearly 
perpendicular to the mean magnetic field fcy <^ k±, and per- 
turbed fields small compared to the mean field SJi <^ Bq. 
The gyrokinetic equation is derived by averaging over the 
particle Larmor motion. This averaging eliminates the fast 
wave and the cyclotron resonances but retains finite Larmor 
radius effects and coUisionless dissipation via the Landau 
resonance. The averaging procedure also eliminates one of 
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the dimensions of velocity in phase space, reducing the di- 
mensionality of the distribution function from 6 to 5, en- 
abling numerical studies of nonlinear gyrokinetic turbulence 
with current computational resources. The validity of the 
gyrokinetic description for the solar wind is explored in de- 
tail in § 3. 

1.2. Turbulence in the Solar Wind 

One of the principal measurements in the study of solar 
wind turbulence is the magnetic field fluctuation frequency 
spectrum derived from tn situ satellite observations. At 
1 AU, the one-dimensional energy spectrum in spacecraft- 
frame frequency typically shows, for low frequencies, a power 
law with a slope of —5/3, suggestive of a Kolmogorov-like 
inertial range [Goldstein et ai, 1995]; a spectral break is typ- 
ically observed at around 0.4 Hz, with a steeper power law 
at higher frequencies, often denoted the dissipation range 
in the literature, with a spectral index that varies from —2 
to —4 [Leamon et al, 1998b; Smith et ai, 2006]. The gen- 
eral consensus is that the —5/3 portion of the spectrum is 
the inertial range of an MHD Alfvenic turbulent cascade, 
but the dynamics responsible for the break and steeper por- 
tion of the spectrum is not well understood. Various expla- 
nations for the location of the break in the spectrum have 
been proposed: that it is coincident with the ion (proton) cy- 
clotron frequency in the plasma [Denskat et ai, 1983; Gold- 
stein et ai, 1994; Leamon et ai, 1998b; Gary, 1999], or that 
the fluctuation length scale (adopting Taylor's hypothesis 
to convert from temporal frequency to spatial wavenumber 
[Taylor, 1938]) has reached either the ion Larmor radius 
[Leamon et al, 1998a, 1999] or the ion inertial length [Lea- 
mon et ai, 2000; Smith et ai, 2001]. The steepening of the 
spectrum at higher wavenumbers has been attributed to ion 
cyclotron damping [Goleman, 1968; Denskat et al, 1983; 
Goldstein et ai, 1994; Leamon et ai, 1998b; Gary, 1999], 
Landau damping of kinetic Alfven waves [Leamon et al., 
1998a, 1999; Leamon et ai, 2000], or the dispersive nature 
of whistler waves [Stawicki et at, 2001; Krishan and Maha- 
jan, 2004; Galtier, 2006]. 

The energy in turbulent fluctuations in the solar wind 
is not distributed isotropically over wavevector space; 
Matthaeus et al. [1990] demonstrated that, at wavenumbers 
k±Pi ~ 10~^, the distribution of energy is anisotropic with 
respect to the direction of the mean magnetic field, creating 
a "Maltese cross" pattern on a contour plot of the mag- 
netic field two-dimensional correlation function. This obser- 
vation has prompted attempts to quantify the amount of en- 
ergy split between nearly parallel modes ( "slab" modes) and 
nearly perpendicular modes ( "quasi- 2-D" modes) [Bieber 
et al., 1996; Leamon et al., 1998b]; results show that 85- 
90% of the power resides in the nearly perpendicular modes. 
Rather than an anisotropic MHD turbulent cascade, Rud- 
erman et al. [1999] propose that phase mixing could be re- 
sponsible for the generation of perpendicular modes. It was 
recently found that the two components of the MHD-scale 
fiuctuations evident in the Maltese cross inhabit separate 
environments; Dasso et al. [2005] decomposed the Maltese 
cross into dominantly perpendicular modes within the slow 
wind and dominantly parallel modes within the fast wind. 
It is thus crucial to distinguish observations in the slow solar 
wind from those in the fast. 

This is not the only difference between the fast and slow 
solar wind that serves as a warning that the dynamics un- 
derlying the turbulence in each system may be distinct. The 
fast and slow components of the wind appear to arise from 
different coronal origins: the steady fast wind originates 
from funnels of open fields lines in coronal holes [Tu et al., 
2005], whereas the unsteady slow wind appears to come from 
temporarily open streamers and from active regions of closed 



magnetic fields [Habbal et ai, 1997; Woo et ai, 2004; Woo 
and Habbal, 2005]. The imbalance between anti-sunward 
and sunward Elsasser spectra can reach nearly two orders of 
magnitude in the fast wind, while the slow wind has a much 
smaller imbalance, from a factor of a few to approximate 
equality [Grappin et al., 1990; Tu et al, 1990]. A detailed 
study of the electron and ion (proton) temperatures in the 
solar wind by Newbury et al. [1998] revealed a nearly con- 
stant mean electron temperature while the mean ion temper- 
ature increased monotonically with the solar wind speed. In 
addition, the sense of the ion temperature anisotropies dif- 
fer, with Till < Ti± in the fast wind and T^n > Ti± in the 
slow [Rasper et al., 2002; Marsch, 2006]. Horbury [1999] 
and Marsch [1999] provide reviews of other salient differ- 
ences between the fast and slow solar wind. 

In this paper, we take heed of the warning by Tu and 
Marsch [1995] of the dangers of averaging measurements 
over a long time period because of the possibility of mixing 
systems with different properties. A variation in parame- 
ters such as the ion plasma beta j3i and the temperature 
ratio Ti/Te can have a strong effect on the kinetic dynam- 
ics of the solar wind, and an imbalance of the sunward and 
anti-sunward wave fluxes may strongly affect the nonlinear 
turbulent transfer of energy, so the fast streams of solar wind 
may evolve quite differently from the slow streams. To pro- 
vide a concrete example of the effect on kinetic dynamics, 
we have calculated the kinetic damping of linear waves with 
k^Pi = 1 using the gyrokinetic^ linear coUisionless disper- 
sion relation [Howes et al, 2006] over the measured range 
of temperature ratios Ti/Te [Newbury et al, 1998] and ion 
plasma beta /3i [Leamon et al, 1998b] (note, however, that 
the Pi range in this paper was not sorted by solar wind 
speed). Figure 2 presents a contour plot of log (7/0;) on the 
{Pi, Ti/Te) plane, where ui is the linear frequency and 7 the 
linear damping rate. Marked on this plot are the parameter 
ranges of the slow (solid shading) and fast (dashed shad- 
ing) solar wind. Immediately evident is that the Landau 
damping"^ is signiflcantly stronger in the slow solar wind. 
It is also clear from Figure 2 that the damping rate varies 
non-trivially with both ion plasma beta /3i and temperature 
ratio Ti/Te. 

1.3. Forward Modeling 

The complexity of Figure 2 underscores the challenge 
of unraveling the underlying physical mechanisms at work 
in the turbulent solar wind from a statistical ensemble of 
measurements taken over a wide range of the key param- 
eters: the solar wind speed Vsw, ion plasma beta /?;, and 
temperature ratio Ti/Te. Unless two of these parameters 
can be held fixed while the third is varied — and impos- 
ing such a constraint would yield a very small number of 
selected measurements — extracting meaningful correlations 
from statistical analyses of data is very difficult. Com- 
pounding this difficulty is the fact that — although current 
and planned missions, such as Cluster, are able to provide 
simultaneous measurements of the solar wind in multiple 
spatial locations — most of our knowledge of the nature of 
turbulence in the solar wind has been gleaned through data 
streams of single-point measurements in time [Fredricks and 
Coroniti, 1976; Bruno and Garbone, 2005]. To interpret the 
data in an attempt to discern the underlying dynamics of the 
solar wind, a wide range of approaches have been employed, 
many of which rest upon reasonable but unproven assump- 
tions. The various approaches to interpreting data range 
from studies of the intrinsic variations of the solar wind pa- 
rameters as well as their evolution with heliocentric distance, 
to detailed examination of the characteristics of fluctuation 
frequency spectra and the anisotropic transfer of energy in 
wavevector space, to analysis of the fluctuations as linear 
eigenmodes of the plasma, to indirect investigation of tur- 
bulent dissipation through thermodynamic measurements. 
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Several recent reviews [Tu and Marsch, 1995; Bruno and 
Carbone, 2005; Marsch, 2006] cover the enormous breadth 
of knowledge that has been gained through these methods 
over decades of study. All of these methods fall under the 
guise of reverse modeling — analyzing observational data in 
an attempt to divine the underlying mechanisms governing 
the evolution of turbulence in the solar wind. 

Even with the wide reaching characterization of the solar 
wind that the above methods have achieved, the mechanisms 
responsible for the ubiquitous break in the magnetic fre- 
quency spectrum between the inertial and dissipation ranges 
or for the varying spectral index in the dissipation range 
have still not been unequivocally identified [Smith et ai, 
2006]. Magnetized turbulence and kinetic plasma physics, 
as well as the starkly contrasting behavior of the fast versus 
the slow wind, contribute to the overwhelming complexity 
of the solar wind, making progress through reverse modeling 
techniques difficult. In light of recent advances in the gen- 
eral understanding of turbulence in magnetized plasmas and 
of the steady gains in computational resources, this paper 
pursues the complementary approach of forward modeling 
[Oughton and Matthaeus, 2005]. 

A forward modeling approach has been fruitfully em- 
ployed to study the spatial transport and spectral evolution 
of large-scale, low-frequency MHD fluctuations in the solar 
wind [Tu et al, 1984; Tu, 1988; Zhou and Matthaeus, 1990a]; 
see the review by Tu and Marsch [1995] for an extensive dis- 
cussion of this well developed field. In addition, the quest to 
understand MHD turbulence has been spearheaded by for- 
ward modeling through detailed numerical simulations [She- 
halin et al., 1983; Matthaeus et al., 1996, 1998; Cho and Vish- 
niac, 2000; Miiller and Biskamp, 2000; Biskamp and Miiller, 
2000; Maron and Goldreich, 2001; Cho et al, 2002; Cho and 
Lazarian, 2003; Vestuto et al., 2003; Oughton et al., 2004; 
Shaikh and Zank, 2006]. To better understand the turbulent 
dynamics of the solar wind at frequencies in the neighbor- 
hood of the observed break in the magnetic energy spectrum, 
a kinetic treatment is essential [Marsch, 2006; Schekochihin 
et al, 2007]. Yet the complexity of this transition regime 
also affords the opportunity to better constrain the under- 
lying physical mechanisms. Forward modeling through nu- 
merical simulation of weakly coUisional turbulent plasmas at 
this transition has just become possible with current com- 
putational resources and codes based on gyrokinetic the- 
ory. The gyrokinetic framework opens up the opportunity 
to study not only the dynamics of the turbulence but also its 
dissipation and the consequent particle heating. We believe 
that, through detailed nonlinear gyrokinetic simulations of 
the turbulent cascade of energy through this regime, the 
mechanisms responsible for the break in the magnetic en- 
ergy spectrum and for the damping of turbulent fluctuations 
and particle heating will at last be identified. 

To facilitate contact between gyrokinetic numerical simu- 
lations, necessarily of limited dynamic range, and the wealth 
of observational data, we construct a turbulent cascade 
model. The model combines the heuristic theories of strong 
turbulence in the MHD Alfven and kinetic Alfven wave 
regimes with the precise predictions of linear wave damp- 
ing. Although this simple model can account for a number 
of observations in the turbulent solar wind, we emphasize 
the importance of nonlinear numerical simulations both in 
justifying the assumptions made by the model and in en- 
abling detailed comparison to measurements that go beyond 
its scope. 

1.4. Outline of the Paper 

In § 2, we first review the theory of critically bal- 
anced, magnetized turbulence both in the MHD and kinetic 
Alfven wave regimes, constructing an analytical model that 



smoothly connects these regimes (§ 2.1). We then proceed 
to use the model to construct magnetic energy spectra con- 
sistent with those observed in the solar wind and to ar- 
gue that the combination of instrumental effects and linear 
damping may be sufficient to explain a wide spread in the 
observed spectral indices in the dissipation range. The valid- 
ity of gyrokinetics in describing turbulence in the balanced, 
slow solar wind, the effective driving scale, the degree of 
anisotropy, evidence of a kinetic Alfven wave cascade from 
satellite measurements and the (un)importance of the ion 
cyclotron resonance are discussed in § 3. The limitations 
of our approach and related previous work are discussed in 
§ 4. Finally, § 5 summarizes our findings, paving the path 
for nonlinear gyrokinetic simulations of turbulence in the 
slow solar wind. 

2. Magnetic Energy Spectra in the Slow 
Solar Wind 

Here we construct a model that follows the nonlin- 
ear transfer of energy from the low driving wavenumber, 
through the inertial range of the critically balanced MHD 
Alfven wave cascade {k±pi <^ 1), connecting to a critically 
balanced kinetic Alfven wave cascade {k±pi 3> 1), and on to 
the highest wavenumbers where the turbulent energy is dis- 
sipated. To extend the cascade into the kinetic Alfven wave 
regime, we choose to follow the magnetic energy rather than 
the kinetic energy of the turbulent fiuctuations because, as 
the ions decouple from the field fiuctuations at scales below 
the ion Larmor radius, their kinetic energy becomes sub- 
dominant. The model assumes that the Alfven wave energy 
fiuxes along the magnetic field in either direction are bal- 
anced; this condition, not generally satisfied in the fast so- 
lar wind, is satisfied in at least some intervals of the slow 
solar wind [Grappin et al., 1990; Tu et al., 1990]. Having set 
up our model in § 2.1, we solve it numerically and present 
the resulting spectra in § 2.2-§ 2.5. In this section, we will 
assume that uo <^ Q.i and that gyrokinetic theory applies. 
This assumption will be investigated in detail in § 3, and 
the limitations of our model are discussed in § 4. 

2.1. Critically Balanced Turbulence and the Turbulent 
Cascade Model 

Consider a homogeneous magnetized plasma permeated 
by a mean magnetic field of magnitude Bo and driven 
isotropically at an outer scale wavenumber fco with veloc- 
ity uq. We assume that both the driving scale and the ion 
mean free path along the magnetic field are much larger than 
the ion Larmor radius, fcoPi -C 1 and Amfpi ^ pi- 

Following Kolmogorov, we adopt the standard assump- 
tion that the flux of turbulent energy (per unit mass and 
volume) at a given scale is entirely determined by the tur- 
bulence at that scale {locality of interactions). Note that 
we make this assumption at all scales. We write the energy 
cascade rate as 

e{ki_)^C-'"^k^Vkbi. (1) 

The magnetic held fluctuation is written in velocity units, 
b\ = SB^{k±)/4:Tmimi, where SB^{k±)/8Tr is the energy 
density of the magnetic held fluctuations perpendicular to 
the mean held integrated over all possible values of fcy . The 
electron fluid velocity perpendicular to the mean magnetic 
fleld is Vk = v±{k±), which is equal to the ion velocity in the 
MHD limit, k±pi ^ 1. All other symbols are defined in Ta- 
ble 1. In principle the dimensionless "Kolmogorov constanf 
Ci is a function of all "local" dimensionless numbers — i.e. 
Ti/Te, Pi, k±pi, k^\/k^, Vk-bk/(ufc6fc), and Lj/(fcx"±), where 
Lo is the local linear frequency. In addition to equation (1), 
we assume that the linear and nonlinear time scales are equal 
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at all scales. This is the critical balance assumption of Gol- 
dreich and Sridhar, 



LU — LUnl{k±) = C2k±Vk- 



(2) 



With this second assumption, the two constants, Ci and C2, 
depend on the slightly reduced set of dimensionless num- 
bers: Ti/Te, Pi, k±pi, k\\/k±, and Vk • bk/(ffebfe). The criti- 
cal balance assumption leads, as we have already explained, 
to the scale dependent anisotropy, fcy ~ fc||(fcx) ^ k±. To 
make further progress we must use the dynamical equations. 
These are reduced MHD (RMHD) in the Alfven wave scale 
range {k±pi < 1) and electron reduced MHD (ERMHD) in 
the kinetic Alfven wave range {k±pi 3> 1) — see Schekochi- 
hin et al. [2007] for a systematic derivation of both sets of 
equations. Linear theory in both scale ranges relates the 
magnetic field and the velocity via 



Vk = ±Q(fcx)bfc, 



(3) 



where the ± refers to the direction of wave propagation along 
B and 



1, fcxpi < 1 

fcxp>/\/A + 2/(1 + r,/r«), fcxp« > 1 



(4) 



We shall assume, in the spirit of the critical balance as- 
sumption, that equation (3) holds in the nonlinear cascade. 
Similarly, the linear frequency in gyrokinetic theory is given 
by [Howes et al, 2006] 



u — ±uj{k±)k^^VA- 



(5) 



where Zj(fcx) = ce{k±) in both asymptotic ranges k±pi <^ 1 
and k±pi S> 1 but not in the transition region k±pi ~ 1. 
The scale invariance of the RMHD and ERMHD equations 
in their respective scale ranges suggests that Ci and C2 are 
independent of k±pi and fcy/fcx within each scale range. 

The turbulent energy is converted to heat via wave- 
particle interactions and collisions. In the linear regime, 
this is manifested as a damping rate, given in the gyroki- 
netic limit by 



7 = 7(A:x)fc||WA. 



(6) 



We will assume that the turbulent heating rate per unit 
mass at each wavenumber k±_ is 276I — in other words, that 
the linear result persists in the nonlinear regime. This al- 
lows us to write a one-dimensional continuity equation for 
the magnetic energy spectrum in perpendicular wavenumber 
space, analogous to the standard simplest spectral evolution 
model used to describe isotropic hydrodynamic turbulence 
[Batchelor, 1953], 



dt 



-k_ 



dk_L 



+ S{kA 



(7) 



where the rate of change of magnetic energy at a given per- 
pendicular wavenumber is given by the three terms on the 
right-hand side: from left to right they are the flux of en- 
ergy in wavenumber space, a source term representing the 
driving of the turbulence (at a low wavenumber fco), and a 
damping term associated with the linear damping rate at a 
given perpendicular wavenumber. In steady state, we write 



k 



dk_L 



= S(fcx) 



or^3/2„ 7(fc±) 

uj[k±) 



(8) 



where we have used equations (1), (2), (5) and (6) to rewrite 
the damping term. Note that because of dissipation the en- 



ergy cascade rate e is not constant but depends on k±. Cal- 
culating ^{k±) /uj{k±) from linear gyrokinetic theory [Howes 
et al., 2006] yields the gyrokinetic cascade model. Equa- 
tion (8) can be integrated if we assume a form for the func- 
tions Ci and C2. Assuming the source S{kj_) is localized at 
ko we have, for k± > ko. 



e(fcx) = eo exp 



-,3/2^ 7(fex) c'fcj. 



where eo is the rate of energy input at fco- 

Using equations (1) and (3), we obtain the one- 
dimensional magnetic energy spectrum 



k ' 



2/3jL-5/3 -2/3 



(10) 



where e is obtained from evaluating equation (9). If dis- 
sipation were neglected, the energy cascade rate would be 
constant e = eo and equations (4) and (10) would lead to 
the familiar relations: EB{kj_) 

-2/3r 



Ciel^^kT^^'-^ in the MHD 
Alfven wave regime and EB{k±_) = Cieg k^''^'^ p~'^^^[Pi + 



2/(1 -I- Te/Ti)]^''^ in the kinetic Alfven wave regime. We 
will assume (as before) isotropic driving at wavenumber fco. 
Thus assuming equations (1), (2), (3), and (5) hold up to 
the driving scale, where A;|| = fcx = fco and kopi ^ 1, we 
obtain eo = Ci^^^C2^kov'^. 

Equations (1), (2), (3), and (5) yield the typical parallel 
wavelength at scale fcx, 



fc||(fcx) 



,2/3,1/30? 



2/3 



1/3 



(11) 



where e is given by equation (9). Note that the normal- 
ized frequency is defined by equation (5) and is calcu- 
lated from gyrokinetics. It is sufficient, given the brutality 
of our assumptions, to make the approximation a(fcx) = 
cJ(fcx) throughout the scale ranges. Again, if dissipation 
were neglected, e = eo and equations (4) and (11) would 
lead to the familiar relations quoted in the introduction: 



1/3,2/3 



in the MHD Alfven wave regime and kn 



kl''-^kY'^p'r'^''-\(5, + 2/(1 + Te/Ti)^/'^ in the kinetic Alfven 
wave regime. 

(Gyro)kinetic theory enters the model we have con- 
structed via the computation of the linear frequency tZ7(fcx) 
(assuming also a = U) in equations (10) and (11) and also 
the linear damping rate 7(fcx) in equation (9). The unde- 
termined Kolmogorov constants Ci and C2 represent a con- 
siderable source of uncertainty, especially because they may 
depend on the plasma parameters Pi and Ti/Te and, around 
k±pi ~ 1, also on fcx. While acknowledging this uncertainty, 
we believe it is useful to look at the implications of our model 
with some fiducial values of the constants. We note that Ci 
simply normalizes the spectrum, so we only need to specify 
the product C\^^C2 in equation (9), which parameterizes 
the overall strength of the coUisionless damping. We use 
the value C\^^C2 = 6, which was originally inspired by the 
values of the constants measured in MHD simulations of 
the Alfven wave regime (Ci ~ 2.5 and C2 — 1.4 [Maron, 
1998]; these values were used also in Quataert and Gruzinov 
[1999]). Without similar numerical guidance in the kinetic 
regime, we take the constants to be the same at all scales. 
While this guess is expedient it is obviously not entirely sat- 
isfactory. Nonetheless, we evaluate equation (9) for e(fcx) 
(numerically) with these assumptions. We will study the 
sensitivity of our results to the values of Ci and C2 in § 2.5. 
The values of these constants will ultimately be determined 
by future gyrokinetic simulations. 

2.2. A Typical Magnetic Energy Spectrum 
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Let us now present steady-state magnetic fluctuation en- 
ergy spectra from the numerical solution of the gyrokinetic 
cascade model constructed in § 2.1. All spectra are normal- 
ized such that the spectrum value is 1 at fc^ = k±i <^ , 
the lowest plotted perpendicular wavenumber. Thus, we 
plot 

EB{k^) = EB{ki_)/EB{k^,) (12) 

where EB{k±i) = {ko/kj_i)'^''^v\/k±iC2- Note that, given 
the assumption uj <^ Qi, the solution to the model of § 2.1 is 
formally independent of the value of the driving wavenumber 
fcopi; this parameter only affects the overall normalization 
of the spectra plotted here. 

In Figure 3 we plot the normalized one-dimensional mag- 
netic fluctuation energy spectrum, EB{k±), as a function of 
wavenumber k±pi for a plasma with ion plasma beta Pi — 1 
and temperature ratio Ti/Te = 1. A plasma of protons and 
electrons is assumed for all the steady-state spectra pre- 
sented in this paper. Two cases are shown: the solution of 
the gyrokinetic cascade model of § 2.1 using the appropri- 
ate gyrokinetic damping rate ^{k±), and, for comparison, 
an undamped case with the damping term artificially set to 
zero, 7(fcx) = 0. The undamped model recovers the spec- 
tral indices of —5/3 in the MHD regime and —7/3 in the 
kinetic Alfven wave regime, as predicted by the arguments 
neglecting dissipation in § 2.1. The damping at k±pi > 1 
in the gyrokinetic model is sufficient to cause the spectrum 
to fall off more steeply than the undamped prediction of 
—7/3. The steady-state gyrokinetic spectrum obtained here 
clearly demonstrates the exponential roll-off characteristic 
of dissipation [Li et ai, 2001]. 

2.3. Effect of Limited Magnetometer Sensitivity 

The magnetic fluctuation spectrum at frequencies above 
the spectral break from in situ solar wind measurements is 
widely interpreted to behave like a power law rather than 
an exponential decay. We suggest here that the power-law 
appearance of the spectrum in this range may be an in- 
strumental effect. The power spectrum of magnetic fluctua- 
tions from satellite measurements is limited by the fiuxgate 
magnetometer sensitivity; this sensitivity limit can, e.g., be 
clearly seen in Figure 6 of Leamon et al. [1998b] at the high 
frequency end of the spectrum. The noise floor of a flux- 
gate magnetometer at frequencies / > 1 Hz is constant in 
units of nT/\/Hz [Bale, 2007]; hence, on a plot of the one- 
dimensional magnetic energy spectrum, this corresponds to 
a constant noise floor at the high frequency end. 

To test the effect of the fiuxgate magnetometer noise floor 
on the appearance of the magnetic fluctuation spectrum 
from our gyrokinetic cascade model, we mock up the instru- 
mental noise by specifying a constant background value to 
the one- dimensional energy spectrum. Inspired by the spec- 
tra in Figure 6 of Leamon et al. [1998b], we choose this noise 
floor to be approximately two to three orders of magnitude 
below the spectrum value at the fcj_pi = 1 break. In Fig- 
ure 4 we plot the spectra for three steady-state solutions of 
our model with the following parameters, chosen to sample 
across the range of solar wind parameters depicted in Fig- 
ure 2: (1) Pi = 0.5, T./Te = 3; (2) Pi = 3, Ti/T, = 0.6; and 
(3) Pi = 0.03, Ti/Te = 0.175. AU cases include the damping 
rate calculated gyrokinetically. Panel (a) shows the spec- 
tra without the added constant sensitivity limit; panel (b) 
adds a constant sensitivity level at two (spectrum 1) or three 
(spectra 2 and 3) orders of magnitude below the spectrum 
value at k_i_pi = 1. The behavior of each spectrum in panel 
(b) in the range k±pi > 1 more closely resembles a power 
law than the exponential roll-off, which would characterize 
the noiseless spectra; the steady-state solutions are well-fit 



by power laws with spectral indices —7/3, —3, and —4. In 
summary, the instrumental sensitivity limit may produce 
spectra that appear to obey a power-law scaling with some 
effective spectral index even though the underlying spectrum 
actually has an exponential roll-off. 

At present there exists little evidence in the literature to 
test this prediction. Search coil magnetometers are more 
sensitive at these high frequencies. Denskat et al. [1983] 
studied the magnetic field power spectra up to 470 Hz us- 
ing data from both fiuxgate and search coil magnetome- 
ters. The search coil data shows power-law behavior up 
to 470 Hz, but the amplitude of the spectral density drops 
about three orders of magnitude in the gap from the highest 
frequency fiuxgate measurement at 2 Hz to the lowest fre- 
quency search coil measurement at 4.7 Hz, suggesting that 
the higher frequency power law is not a continuation of the 
lower frequency spectrum. Alternatively, the mismatch may 
be due to difficulty performing a relative calibration of the 
instruments, but without an overlap in frequency the dif- 
ference cannot be resolved. One recent paper [Alexandrova 
et al., 2008] shows search coil measurements of magnetic 
field spectra up to 10 Hz, a factor of 20 in frequency above 
the breakpoint; for the 57-minute interval of the fast solar 
wind studied, the dissipation range spectrum shows neither 
an exponential roll-off nor evidence of instrumental noise 
up to 10 Hz. During one 7-minute interval, the dissipation 
range index is —2.33; during another 7-minute interval, it 
is —2.50; and over the whole 57-minute interval, it is —2.6. 
These results are not at odds with our model predictions 
for weak damping — the dissipation range spectrum is ex- 
pected to appear as a —7/3 power law to higher frequen- 
cies before the damping becomes strong enough to cause 
the spectrum to roll off exponentially. We must also keep in 
mind that the fast wind generally does not have a balance of 
Elsasser energy fluxes, an assumption of the present simple 
cascade model. An interesting aspect of the measurements 
in Alexandrova et al. [2008] not discussed in their paper is 
the time variability of the spectral index in the dissipation 
range. There is a clear need for a detailed study of solar 
wind turbulence at frequencies above 1 Hz, taking into ac- 
count the varying plasma parameters, the imbalance in the 
Elsasser energy fluxes, the sensitivity limits of the various 
instruments, and the time variability of the dissipation range 
spectral index. 

It is instructive to consider what we will learn if such a 
study does not flnd magnetic energy spectra demonstrating 
a slow, exponential roll-off. The absence of an exponential 
decay could indicate that using the linear damping rates in a 
turbulent plasma is incorrect. It is important to emphasize 
that the predictions of gyrokinetic theory for solar wind tur- 
bulence are not equivalent to the predictions of this simple 
cascade model — nonlinear gyrokinetic simulations of turbu- 
lence will play a key role in testing this model, assessing the 
robustness of the predicted exponential fall-off, and refining 
the model as necessary. 

2.4. The Effective Spectral Index in the Dissipation 
Range 

Figure 4 shows that, over the range of the plasma param- 
eters Pi and Ti/Te measured in the solar wind, the observed 
wide spread of the spectral index in the dissipation range 
from —2 to —4 [Leamon et al., 1998b; Smith et al., 2006] 
can be explained by the damping of kinetic Alfven waves via 
the Landau resonance. If the Landau damping is negligible, 
the spectral index is expected to be —7/3, close to the ob- 
served upper limit; over the range of parameters relevant to 
the solar wind, our cascade model gives a lower limit to the 
(effective) spectral index of about —4, consistent with obser- 
vations. The measured effective spectral index for fcj_pi > 1 
is a complicated function of the parameters Pi, T/Te, and 
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kopi as well as of the sensitivity limit of the magnetome- 
ter. If it is indeed the linear damping that explains the wide 
spread in the measured values, the complexity of this pa- 
rameter dependence is a natural outcome of the nontrivial 
behavior of the coUisionless damping as a function of Pi and 
Ti /Te over the range of parameters in the slow and fast solar 
wind depicted in Figure 2. Plots of the measured spectral in- 
dex versus a single plasma parameter are unlikely to reveal 
strong correlations due to variation in the other parame- 
ters within the data set. To make headway here, a forward 
modeling approach may prove more fruitful. One needs to 
incorporate all of the measured parameters to determine the 
predicted spectral index for an interval of solar wind data. 

2.5. Uncertainty in the Kolmogorov Constants 

The spectrum derived from our cascade model depends 
on the values of the Kolmogorov constants in the theory. 
As we noted above, while Ci merely normalizes the spec- 

3/2 

trum, the product Cj^ C2 affects the competition of the 
nonlinear transfer of energy with the linear damping [see 
equation (9)]. If this product is increased (decreased), this 
effectively makes the linear dissipation stronger (weaker) rel- 
ative to the nonlinear transfer. These order unity constants 
must be determined from nonlinear simulations or experi- 
mental data. To judge the effect of varying these constants 
on the steady-state solution. Figure 6 presents the gyroki- 
netic spectrum for Pi = 1 and Ti/T^ — 1, with this product 
taken to be Cl^^C2 ~ 3,6, and 12. The resulting variation 
in the effective spectral index is -1-1/3 or —2/3, so all of our 

steady-state spectra suffer this systematic uncertainty. In 

3/2 

the rest of this paper we use the fiducial value C-^ C2 = 6. 

3. Applicability of Gyrokinetics to the Slow 
Solar Wind 

In this section, we argue that gyrokinetics is applicable to 
the dynamics of the turbulence in the balanced, slow solar 
wind at the dissipation scale. Gyrokinetics is not, of course, 
applicable at the driving scale L — 27r/fco, where the mag- 
netic field fluctuations are of the order of the mean field and 
are relatively isotropic. Formally, as described in § 1.1, a 
number of conditions must all be met for the gyrokinetic 
approach to be valid: pi <ti L, fcy <C k±, SS -ti -Bo, and 

^ Sli. Let us consider whether each of these conditions 
is satisfied in the solar wind. In § 3.1 we estimate the effec- 
tive driving scale {pi <^ L). An estimate of the anisotropy 
k\\ ^ k± at the dissipation scale k±pi ~ 1 follows in § 3.2. 
The observed amplitude of the fluctuations 5B <^ Bq is 
discussed in § 3.3. In § 3.4 we present wave phase velocity 
evidence from satellite measurements that is consistent with 
a kinetic Alfven wave cascade but inconsistent with the on- 
set of ion cyclotron damping — indirect evidence for uj <^ Q,i. 
In § 3.5 we show (in some detail) that the turbulent frequen- 
cies are usually less than the cyclotron frequency {uj <^ Qi). 
Finally, in § 3.6, we explain what happens if the cyclotron 
resonance is reached in the dissipation range. 

3.1. Turbulence Driving Scale 

The nature of the turbulent fluctuations in the dissipa- 
tion range of the solar wind depends on the driving scale L 
at which energy is injected into the turbulent cascade. The 
anisotropy predicted by the critically balanced Alfven wave 
cascade in GS theory is scale dependent, so this raises two 
questions: is the solar wind driven isotropically or anisotrop- 
ically, and what is the effective driving scale L? Unfortu- 
nately, the mechanisms responsible for driving turbulence in 
the solar wind are unknown; the matter is further compli- 
cated by the radial expansion of the solar wind as it travels 
outward from the sun and the measured imbalance of anti- 
sunward vs. sunward directed Alfven wave flux [Grappin 



et ai, 1990; Tu et ai, 1990; Tu and Marsch, 1990]. Never- 
theless, we can make some headway on this complex problem 
if we ignore the speciflc details of the driving mechanism and 
the radial expansion and simply estimate the anisotropy of 
the balanced, slow solar wind turbulence from in situ mea- 
surements at 1 AU. 

As discussed in § 1.1, determining the isotropic driving 
wavenumber kopi — 2tt/L fixes the path of the critically 
balanced cascade as depicted in Figure 1. Equivalently, a 
measurement of the wavenumber anisotropy at any scale 
in the solar wind can also be used to fix the critical bal- 
ance relation; here we use this measurement to estimate 
the resulting effective value of kopi. From measurements of 
the anisotropy of the spatial correlation of magnetic field 
fiuctuations in the slow solar wind, presented in the left 
panel of Figure 1 from Dasso et al. [2005], we estimate 
the anisotropy fcx/fcy = r^^/r± ~ 1.5 at an average scale 
2-K/k± = 6 X 10^*^ cm. If we assume that the turbulence 
is strong at the isotropic driving scale and that the energy 
cascade rate is a constant e = eo between the driving and 
the observed scales (a reasonable assumption given that ob- 
served spectra consistently show a spectral index of —5/3 
[Smith et ai, 2006] and that the damping rate for Alfven 
waves is small, 7 ^ a;, at these scales), we can estimate 
the driving scale from which the observed turbulence arose 
using equation (11), 

kopi = (jT-^ k±pi. (13) 

which applies for scales in the MHD Alfven wave regime 
k±pi <^ 1. Using a mean ion (proton) temperature for the 
slow wind Ti ~ 4.5 x 10* K [Newbury et at, 1998] and the 
geometric mean of the range of magnetic field magnitude in 
Leamon et al. [1998b] of -Bo — 10~* G, we obtain the ion 
Larmor radius pi ~ 3 x 10® cm. Hence, we find kopi ^ 10~*, 
or an isotropic driving scale -L ~ 2 x 10^^ cm. 

This value is consistent with the estimate by Matthaeus 
et al. [2005] of the mean magnetic field correlation length 
in the solar wind (over all solar wind speeds) of L ~ 
1.2 X 10^^ cm. Figure 3 of Tu and Marsch [1990] presents 
spectra of the outward directed Elsasser flux of Alfven waves 
measured by the Helios spacecraft missions which show, for 
the slow solar wind at 0.9 AU,* the beginning of the —5/3 
spectrum at a reduced wavenumber fc/27r ~ 10^^ km^^, cor- 
responding to a driving scale of -L = 10^^ cm. Finally, an es- 
timate of the correlation length of the anti-sunward Elsasser 
held at 0.3 AU is L = 1.6 x lO" cm [Tu and Marsch, 1995]; 
scaling this length for radial expansion to 1 AU gives a value 
of L = 5 X 10^^ cm. These four independent determinations 
of the driving scale, at which we assume the fluctuations 
are isotropic, are consistent to within an order of magni- 
tude. This suggests that the driving is at least approxi- 
mately isotropic and that kopi = 10"* is a good choice for 
a flducial value of the isotropic driving wavenumber. 

It is perhaps appropriate to discuss why this number 
is physically sensible. The simplest interpretation of so- 
lar wind turbulence is that energy resides in large-scale, 
energy containing modes driven in the solar corona, and 
those modes are subsequently left to decay freely in the wind 
within the inner heliosphere. A reasonable estimate for the 
decay time of turbulent energy characterized by a velocity 
v±_ at scale fcx is the eddy turnover time r ~ {k_i_v±_)~^ . For 
the driving scale estimated above, L = 27r/fcx = 2x 10^^ cm, 
and a turbulent velocity estimated from the slow solar wind 
Elsasser spectra in Figure 4 of Grappin et al. [1990] of 
= 2 km/s, this decay time is r ~ 50 h. The travel 
time for slow solar wind with a velocity Vsw ~ 350 km/s to 
reach R = 1 AU is t = R/vsw ~ 100 h. If there is energy 
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in much larger scales than our estimated value for L, those 
larger scale fluctuations will not have had time to transfer 
their energy in the nonlinear cascade. On the other hand, 
energy in modes much smaller than L in the solar corona 
would have already cascaded to the dissipation scale and 
would not be observed at 1 AU. Indeed the near agreement 
between the estimated turbulent decay time, r ~ 50 h, and 
the travel time of the slow solar wind to 1 AU, t ~ 100 h, 
is probably not coincidental. It suggests that, within 1 AU, 
the turbulence in the slow solar wind is freely decaying and 
that we are seeing the energy injected at the largest scale 
that can cascade on the travel time. 

3.2. Anisotropy at the Ion Larmor Radius Scale 

Assuming an ion number density in the slow solar wind 
of rii ~ 20 cm~'^ Leamon et al. [1998b] and an electron 
temperature of Te ~ 1.3 x lO'^ K [Newbury et al, 1998], 
we find the ion plasma beta (3i ~ 0.3 and temperature ra- 
tio Ti/Te — 0.35; for these parameters, the normalized fre- 
quency from the gyrokinetic dispersion relation at fcj_pi = 1 
is a; ~ 1.45. Substituting these values, along with our fidu- 
cial isotropic driving wavenumber kopi = 10~*, into equa- 
tion (11) and taking for simplicity e — eo, we obtain the 
estimate k\\/k± ~ 0.04 at k±pi = 1. As will be demon- 
strated in Figure 8, by neglecting dissipation, this estimate 
represents an upper limit on the value of fcy/fcx. 

3.3. Fluctuation Amplitude in the Dissipation Range 

One might be concerned that gyrokinetics requires a 
strong mean magnetic field in order for the gyroaveraging 
procedure to be valid, and that the solar wind does not con- 
tain such a strong mean field. At the driving scale L of the 
turbulence in the solar wind, the magnetic field has fluctu- 
ations of the same order as the mean field, 5B ~ Bq. How- 
ever, as was realized by Kraichnan [1965] , even in a plasma 
with no mean magnetic field, the magnetic field of the large- 
scale fluctuations can serve as an effective mean field for 
smaller-amplitude fluctuations on smaller scales. Raw mag- 
netometer data, for instance Figure 3-5a of Tu and Marsch 
[1995] or Figure 1 of Bale et al. [2005], show order one fluctu- 
ations of the magnetic field vector components, correspond- 
ing to large changes in the direction of the mean magnetic 
field (the magnitude of the field has much smaller fractional 
fluctuations) on timescales from 15 minutes to 1 hour. On 
shorter timescales, however, the field direction remains rel- 
atively constant with small amplitude, short-timescale fluc- 
tuations about the mean. This is in agreement with Kraich- 
nan's hypothesis, showing fluctuations of both small ampli- 
tude and scale about a mean magnetic field arising from the 
larger scale fluctuations — these small fluctuations are pre- 
cisely those which are well-described by gyrokinetics. Once 
the cascade has progressed to scales a couple of orders of 
magnitude smaller than the driving scale, the amplitude 
of the fluctuations becomes small compared to the local 
mean field amplitude and the fluctuations become signifi- 
cantly anisotropic, 5B/_Bo <gi 1 and /cy <gi k±. Gyrokinetics 
then becomes a valid description of the fluctuations, and the 
deeper one goes into the inertial range, the better it works. 
From equation (10), taking e = eo, we have 

^ ^ ^ ^ / fcc, (14) 

Bo VA C2\k^J • ^ ' 

Thus, for our estimated ko we have S'B/Bo < 0.03 for 

k±pi 1- 

3.4. Evidence against the Cyclotron Resonance 

There exists evidence from in situ satellite measurements 
consistent with the existence of a kinetic Alfven wave cas- 



cade, at least in some cases. Wave phase velocity informa- 
tion in the region of the spectral break can be used to dis- 
tinguish the characteristics of the underlying wave modes. 
Bale et al. [2005] have produced the first simultaneous mea- 
surements of the fluctuating electric and magnetic fields, 
demonstrating that the wave phase velocity above the spec- 
tral breakpoint increases, as shown in panel (b) of Figure 3 
in Bale et al. [2005]. Using the Taylor hypothesis to con- 
vert the frequency spectrum to a wavenumber spectrum, 
the spectral breakpoint occurs near the ion Larmor radius, 
kpi ~ 1. 

Using the full Vlasov-Maxwell linear dispersion relation 
for a hot collisionless plasma (see Quataert [1998] for a de- 
scription of the code that we use to solve it), we plot in 
Figure 7 the normalized wave phase velocity uj/k^^VA vs. the 
normalized wavenumber kpi. Plasma parameters are cho- 
sen to be /3i = 1 and Ti/Te = 1. Although in the gyroki- 
netic theory ui/k^^VA = Lj{k±) is a function of k± alone, in 
hot plasma theory cu/k^^VA depends also on fcy. The par- 
allel wavenumber needed for this calculation is determined 
by using equation (11) with no dissipation, e = eo. The 
change in wave phase velocity is plotted for driving scales 
kopi = 10~®, 10~*, 10~^, and 1. The wave phase veloc- 
ity increases until the ion cyclotron frequency is reached; 
then, the frequency asymptotes uj —> Qi as ion cyclotron 
damping becomes important (see Appendix A), so uj/k\\VA 
decreases due to the increasing fc||(fcx) [see equation (11)]. 
The measured wave phase velocity in Figure 3 of Bale et al. 
[2005] does not show this downturn, suggesting that the ion 
cyclotron damping remains unimportant below the highest 
wavenumber measured, k^pi ~ 10. In the limit fcy ^ fcx, 
the upturn of the wave phase velocity above the breakpoint 
in the magnetic fluctuation spectrum is consistent with un- 
derlying wave modes consisting of kinetic Alfven waves with 
uj <^ ^i. We do note, however, that this increase of phase 
velocity may also be consistent with fast MHD/whistler 
mode of Vlasov-Maxwell kinetic theory for certain wavevec- 
tor anisotropies. For either case, though, it appears that the 
break is not due to the onset of ion cyclotron damping. 

3.5. Conditions for Reaching the Ion Cyclotron 
Resonance 

Although the Bale et al. [2005] measurements of the elec- 
tric fields in the solar wind are consistent with kinetic Alfven 
waves with frequencies lo <^ Q.i, there is also indirect evi- 
dence that fluctuations at or above the cyclotron frequency 
are present in the solar corona and solar wind. For example, 
minor ions such as O"""^ are selectively heated to high per- 
pendicular temperatures [Kohl et al., 1997, 1998]. Also, the 
adiabatic invariant p, oc T^/B for ions and alpha particles 
is observed to increase with heliocentric distance [Marsch 
et al., 1983], indicating some mechanism for perpendicu- 
lar heating of these species, particularly in the fast wind. 
Since (the collisionless part of) gyrokinetics preserves the 
adiabatic invariant p, it cannot describe such heating. This 
heating may come from fast wave fluctuations [Chandran, 
2006] or velocity-space instabilities [Scarf et al., 1967; Evi- 
atar and Schulz, 1970; Gary et al., 1976; Kasper et al, 2002; 
HelUnger et al., 2006; Marsch, 2006] rather than the Alfven 
wave cascade, but it is nevertheless important to know where 
in (fex,fc||) space the turbulent cascade may start to be af- 
fected by the cyclotron resonance. In this section we employ 
the steady-state solutions of the turbulent cascade model 
presented in § 2 to estimate the perpendicular wavenumber 
at which the cascade reaches the cyclotron frequency — this 
is the same wavenumber at which a gyrokinetic treatment 
fails to be valid. 

From equation (5) and equation (11) we obtain 

^ = ^aJ(fcx) = [k^p.r\kop.fr^^ (±] (15) 
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where e/eo is given by equation (9) and we have assumed 
a = tiJ(fcx) for all A:^. In Figure 8, we plot uj/Q.i and fcypi 
for both the undamped and gyrokinetic cascade models with 
Pi — 1, TijT^ = 1, and fcoPi = 10~*. It is unknown if the 
linear damping rates are a good estimate of the damping of 
nonlinear turbulent fluctuations, but the undamped model 
certainly provides upper limits for the frequency and paral- 
lel wavenumber. Figure 8 shows that the undamped model 
reaches the cyclotron frequency at fcxPi — 12, whereas 
the model using the linear gyrokinetic damping rates never 
reaches the cyclotron frequency. The plot of k\\pi vs. k±pi 
for the gyrokinetic model in Figure 8 reaches a maximum at 
fcxpi ~ 5 and then fcypi begins decreasing. In this typical 
case, it is clear that, if the cyclotron frequency is reached, it 
is only after the fluctuations have already cascaded to a high 
perpendicular wavenumber, k±pi ^ 1, at which point most 
of the energy in the cascade has already been dissipated. 
Consequently, the gyrokinetic treatment will be valid over a 
large part, if not all, of the dissipation range. 

The decrease in fcypi seen in the solution of the gyroki- 
netic model in Figure 8 is likely an unphysical result arising 
from the imposition of critical balance, a condition requiring 
that the turbulence remains strong. The kinetic damping of 
the turbulent fluctuations becomes significant at this point, 
so the nonlinear frequency begins to fall below the linear 
frequency. Instead of fcypi decreasing, it is more likely that 
the turbulence becomes weak, a condition not treated within 
our simple cascade model. The realm of weak turbulence, 
LJ > uj„i, exists above the line of critical balance as depicted 
in Figure 1. In weak turbulence theory, the cascade to higher 
parallel wavenumber is suppressed [Goldreich and Sridhar, 
1997; Ng and Bhattacharjee, 1997; Galtier et al., 2000; Lith- 
wick and Goldreich, 2003], so a more physical result is that 
the value of fcypi should remain constant, rather than de- 
crease, as the cascade progresses to higher k±pi. Adjusting 
the model to prevent the value of fcypi from decreasing re- 
sults in a sharper cutoff in the magnetic energy spectrum 
at the high wavenumbers where the parallel wavenumber is 
held constant — this is the result of a decrease in the ratio of 
the nonlinear transfer frequency to the linear damping rate 
as the turbulence weakens. The simplicity of our cascade 
model limits our ability to determine the correct physical 
behavior in this regime; nonlinear simulations are necessary 
to determine accurately the behavior of the turbulent cas- 
cade as the kinetic damping becomes significant. 

Analogous to our gyrokinetic cascade model, we can con- 
struct a hot plasma cascade model using the frequencies and 
damping rates from the linear hot plasma dispersion rela- 
tion [Stix, 1992; Quataert, 1998] instead of the gyrokinetic 
values. In hot plasma theory, however, the normalized wave 
frequency u/{k\'^VA) is no longer a function of k± alone, but 
depends also on fcy . Since the linear hot plasma and gy- 
rokinetic theories agree when the gyrokinetic approximation 
is valid, we will continue to use the solution of the gyroki- 
netic model to specify the path of the cascade in the (fcx , fey ) 
plane. This is probably at least qualitatively correct because 
cyclotron damping does not in fact become large. 

To quantify the conditions where the cyclotron resonance 
terms becomes important we define a threshold value of k± 
to be the value along the cascade path at which the solu- 
tions of the gyrokinetic and hot plasma dispersion relations 
diverge by 10%. This threshold is plotted (solid line) in 
Figure 9 for a range of fii with parameters Ti/T^ = 1 and 
kopi = 10~*. We probe the sensitivity of the threshold k±_ 
to our assumptions by plotting the same 10% threshold for 
the path in the (fcx,fc||) plane calculated from the gyroki- 
netic model with no damping (dashed line) . The dotted line 
is an estimate of this threshold determined by solving the 
fcxPi ^ 1 limit of equation (15) with ui/Jli = 1 and e = eo 
to find 



LPi 



[kopi) Pi 



1 + 



Pi{i + n/T,) 



1/4 



(16) 



(this rough estimate is only plotted over the typical range 
of Pi in the solar wind near Earth). 

The solid line in Figure 9 represents our best estimate 
of the threshold at which the cyclotron resonance becomes 
important and the gyrokinetic approximation fails. If the 
effective damping of the fluctuations in the nonlinear turbu- 
lence is less than the linear damping rates, the true thresh- 
old will fall within the vertically shaded region. This region 
is bounded from below by the short-dashed line — our esti- 
mate of the threshold with no damping — and represents the 
most conservative estimate of the limit of applicability of gy- 
rokinetics. Thus, throughout the diagonally shaded region, 
gyrokinetics is certainly a valid description of the turbulent 
dynamics in the plasma, and this region covers a substan- 
tial part of the parameter range of interest in the study of 
turbulence in the solar wind. 

The cyclotron resonance has a strong effect on the cas- 
cade only when it is reached at a relatively low wavenumber, 
k±pi < 1; otherwise, the cascade is likely to be rather heav- 
ily damped before reaching the cyclotron frequency, leav- 
ing little energy available for perpendicular heating. Equa- 
tion (16) provides a rough estimate of how the wavenum- 
ber threshold in the undamped case scales with the plasma 
and system parameters. This undamped threshold has lit- 
tle dependence on the temperature ratio Ti/Te- Thus, the 
cyclotron resonance occurs early on in the dissipation range 
only in plasmas with Pi < 0.01 or when the isotropic driv- 
ing wavenumber is much larger than our fiducial estimate 
koPi > 10"'' . 

3.6. Effect of the Ion Cyclotron Resonance 

Let us now discuss what happens if the cyclotron res- 
onance is reached in the dissipation range {k±pi > 1) and 
before the kinetic Alfven wave cascade is converted into heat 
by Landau damping (fcxPe < 1). The idealized asymptotic 
version of the appropriate linear theory (the case k±pi ^ 1, 
k±ps ^ 1 at the cyclotron resonance) can be worked out an- 
alytically from the general hot plasma dispersion relation — 
this is done in Appendix A. The key points that this calcu- 
lation demonstrates are 

• the cyclotron effects only matter when the frequency of 
the fiuctuations is very close to one of the resonances: for 
integer n> 1, {uj- uQi)/ (\kii\vAk±pi) ~ l/^/k±p^■, 

• away from the resonances, the ions effectively have 
a Boltzmann response and kinetic Alfven waves (weakly 
damped by the electron Landau resonance) exist both be- 
low the ion cyclotron frequency {uj < fli) and between the 
resonances {nQi < u < {n + 1)0^); 

• at the resonances, there is a conversion between kinetic 
Alfven waves and ion Bernstein waves — electrostatic waves 
whose frequencies are close to nili [cf. Li and Habbal, 2001]; 

• because k±pi ^ 1, the ion cyclotron damping is expo- 
nentially weak (~ e"'"^''') and very localized in frequency. 

These points are illustrated in Figure 10. Together, they 
imply that the effect of passing through the cyclotron reso- 
nance on the kinetic Alfven wave cascade is unlikely to be 
drastic: in order to be subjected to even a weak cyclotron 
damping, the ffuctuations have to couple nonlinearly into a 
rather narrow frequency band around fii; it appears more 
likely that most of the turbulent energy would leapfrog the 
resonance altogether and continue as a kinetic Alfven wave 
cascade, until converted into heat by electron Landau damp- 
ing (at fcxPe ~ 1). However, it is probably prudent to ac- 
knowledge that linear calculations such as the one presented 
in Appendix A may have limited utility in predicting non- 
linear behavior. 
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It is important to understand that the calculation in Ap- 
pendix A considers an idealized asymptotic solution. It 
is known that the precise mode conversion and damping 
properties at the cyclotron resonances depend rather sensi- 
tively on the plasma parameters: such behavior at oblique 
(fcj_ > fc||) propagation was explored in detail by Li and 
Habbal [2001]. To demonstrate some of the non-asymptotic 
possibilities that may be present in the solar wind, it is per- 
haps useful to consider here a rather extreme low-/3i example 
in which the ion cyclotron resonance is reached quite close 
to the beginning of the dissipation range. 

In Figure 11, we present the cascade path on the (fcx, fcy) 
plane (solid) and the frequencies along that path calculated 
by numerically solving both the linear gyrokinetic (dashed) 
and hot plasma (dotted) dispersion relations for plasma pa- 
rameters Pi — 0.03, Ti/T^ — 3, and fcoPi = 2 x 10"^ (chosen 
to highlight mode conversion features) . Two branches of the 
hot plasma dispersion relation are plotted: in the MHD limit 
k^Pi <^ 1, these correspond to the MHD Alfven wave and 
MHD fast wave. As discussed in § 1.1, the MHD Alfven wave 
transitions to a kinetic Alfven wave at fcxPi > 1- The cor- 
responding solution of the gyrokinetic dispersion relation, 
which does not include the cyclotron resonance, continues 
right through = fi^. The effect of the cyclotron resonance 
on the Alfven branch of the hot plasma solution is that, 
as Q,i, the kinetic Alfven wave converts into an ion 

Bernstein wave at k±pi ~ 3. The fast wave has a similarly 
interesting behavior in this example. At k±pi ~ 0.3, the 
fast wave converts to an ion Bernstein wave; eventually this 
same branch, at k±pi ~ 5 connects onto a continuation of 
the kinetic Alfven wave in the range Qi < u> < 2^1i before 
becoming an n = 2 ion Bernstein wave at fcxPi — 8. 

In a turbulent cascade, different wave modes can couple 
nonlinearly, but if the nonlinear transfer is local in wavenum- 
ber, the coupling will be strong only when those modes have 
similar wavevectors and frequencies. Thus, in our exam- 
ple, at k±pi ~ 0.1, the two branches (the fast and Alfven 
modes) have frequencies that differ by an order of magni- 
tude, so the coupling is likely to be weak; but, at k±pi ~ 4, 
their frequencies differ by only a factor of 2, so the coupling 
may be much more efficient. In this case, the coupling may 
enable the turbulent energy of the kinetic Alfven wave on 
the low-frequency branch to cascade nonlinearly to the con- 
tinuation of the kinetic Alfven wave on the high-frequency 
branch. The efficiency of such a process and the amount of 
perpendicular cyclotron heating that can occur as the cas- 
cade passes through the resonance is likely to depend on the 
damping rates of the modes involved. In Figure 12 we plot 
the damping rates onto ions and electrons for the solutions 
of the gyrokinetic (dashed) and hot plasma (dotted) dis- 
persion relations for the same parameters as in Figure 11. 
In the gyrokinetic case, the ion Landau damping falls be- 
low the plotted range, while the electron Landau damping 
is significant (indeed, it is this damping that served as the 
primary kinetic input in our cascade model). For the low- 
frequency (Alfven ) branch (upper panel), there is a rise of 
ion cyclotron damping at fcxPi — 3 and a slight decrease in 
the electron Landau damping. For the high-frequency (fast) 
branch (lower panel), there is a narrow peak of ion cyclotron 
damping at k±pi ~ 0.16, the point where the fast wave fre- 
quency crosses uj — Hi. Two more ion cyclotron damping 
peaks occur in wavenumber ranges where the fast branch 
takes on the characteristics an ion Bernstein wave. (Similar 
cyclotron damping peaks for the low- and high-frequency 
branches are seen in the more asymptotic case considered in 
Appendix A; see Figure 10). Thus, if the kinetic Alfven wave 
cascade reaches the ion cyclotron frequency, it is possible for 
some of the energy to go into perpendicular ion heating via 
coupling to the ion Bernstein wave modes. However, for 
parameters typical of the slow solar wind, and even for the 
rather extreme low-/3i example considered here, the ion cy- 
clotron damping is weak because it only sets in at k±pi > 1; 



in Figure 12 it remains at least an order of magnitude below 
the electron Landau damping. 

As noted above, when the cascade reaches the ion cy- 
clotron frequency, the path of the cascade on the (fcx,fc||) 
plane, calculated using the gyrokinetic cascade model, is 
unlikely to be quantitatively correct. Nonetheless, we can 
tentatively examine the differences in the resulting magnetic 
energy spectra between the gyrokinetic and hot plasma cas- 
cade models (defined in § 3.5). These are plotted in Fig- 
ure 13 for the same parameters as the example in Figures 11 
and 12. Because the ion cyclotron resonance does not signif- 
icantly alter the total damping rates (see the upper panel of 
Figure 12), the resulting spectra difi'er very little, both pro- 
ducing an effective dissipation range spectral index near —4. 
This suggests that the effect of the ion cyclotron resonance 
cannot be identified based solely on the magnetic fluctuation 
energy spectrum; an examination of the underlying dynam- 
ics in concert with a forward modeling approach, like the 
wave phase velocity study in § 3.4, is necessary if one hopes 
to distinguish between ion cyclotron and electron Landau 
damping and, consequently, between ion and electron heat- 
ing. 

4. Discussion 

The gyrokinetic cascade model developed in this paper is 
meant to describe in a simple way the physical mechanism — 
the nonlinear cascade through the MHD Alfven wave and 
kinetic Alfven wave regimes — by which the vast reservoir of 
energy contained in large-scale turbulent eddies in the solar 
wind cascades to high wavenumber and frequency and ulti- 
mately dissipates due to coUisionless damping. Our aim is 
to reconcile observational data with theoretical ideas about 
the energy transfer to the small dissipative scales, and even- 
tually to thermal energy of the plasma particles. In that 
light, this section contains a lengthy discussion of the limi- 
tations of our model, followed by a discussion of its relation 
to previous work. 

4.1. Limitations of Our Model 

One of the key assumptions of our model is that the tur- 
bulence can be treated as if it were driven isotropically at 
some large scale, characterized by an effective isotropic driv- 
ing wavenumber kopi- Two exceptions can be made to this 
assumption: the turbulence may not be driven isotropically, 
and it may be driven at more than one scale. 

Possibility of Anisotropic Driving: To dispense with 
this question, we recall that, as we noted in § 1.1, our 
model does not strictly require that the turbulence is driven 
isotropically. The critical balance path in the {k±, fcy) plane, 
as plotted in Figure 1, simply requires a single point to fix 
the relation — we use the parameter koPi as a convenient 
characterization of that point. The anisotropy of the slow 
solar wind turbulence at k±pi ~ 10~^, measured by Dasso 
et al. [2005], supplies the necessary data point to pin down 
the scale-dependent anisotropy of the GS theory; we de- 
termine the characteristic value of of fcopi based on that 
point. A strong test of this model would be a measurement 
of the anisotropy, similar to that in Matthaeus et al. [1990], 
at a smaller scale around fcx Pi ~ 0.1. At this wavenum- 
ber the anisotropy is predicted to be much stronger with 
fcx/fc|| ~ 10. 

A note is, perhaps, necessary regarding our use of fre- 
quencies and damping rates from the gyrokinetic dispersion 
relation for the entire turbulent cascade even though gyroki- 
netics is obviously not valid at the isotropic driving scale. 
While our model undoubtedly does not describe the driving 
scale well, the effect of driving is probably well modeled by 
the simple idea of a constant energy input: the damping rate 



HOWES ET AL.: DISSIPATION OF THE SOLAR WIND 



X - 11 



is negligible at the driving scale and the Alfven frequency 
gives a good order-of-magnitude approximation to the non- 
linear cascade rate, so the values from the gyrokinetic theory 
are satisfactory. As we noted in § 3, the assumptions behind 
the gyrokinetic theory are ever better satisfied as we proceed 
deeper into the inertial range. 

Possibility of Small-Scale Driving: The second ex- 
ception is that energy may be injected at multiple scales. In 
§ 3.1, we argued that our estimated effective driving scale 
was consistent with the turbulence driven at the Sun and 
then decaying in transit to 1 AU. It is unlikely that there 
are large-scale energy sources affecting the turbulence dur- 
ing this transit. It has, however, been suggested that Kelvin- 
Helmholtz instabilities arising from shear between fast and 
slow solar wind streams could continually drive turbulence in 
the solar wind [Coleman, 1968; Roberts et al, 1987; Roberts 
et al., 1991; Roberts et al, 1992; Zank et al, 1996]; such a 
mechanism could inject energy at scales smaller than our 
estimated isotropic driving scale. However, magnetic fluc- 
tuation frequency spectra for the slow solar wind typically 
show a very clean —5/3 spectrum all the way to the low 
frequencies associated with our fiducial driving wavenumber 
kopi ~ 10^" [Tu and Marsch, 1990; Grappin et al., 1990; Tu 
and Marsch, 1995], suggesting a continuous inertial range 
below the energy injection scale. If significant energy were 
injected at smaller scales, one would expect to see structure 
in the magnetic energy spectrum at the associated frequen- 
cies. The lack of such structure argues against significant 
energy injection at smaller scales. Hence, we conclude that 
a single effective isotropic driving wavenumber koPi is a rea- 
sonable assumption. 

Plasma Instabilities: Another potential source of 
small-scale energy injection neglected in our model is the 
possible presence of temperature anisotropy instabilities, 
e.g., the mirror or firehose instabilities. Adiabatic expan- 
sion of the solar wind drives ions towards Ti± < T^n, which 
should trigger firehose instabilities [Scarf et al, 1967]; in- 
deed, there exists observational evidence that the solar wind 
ion temperature anisotropies are constrained by the fire- 
hose marginal stability criterion [Eviatar and Schulz, 1970; 
Kasper et al, 2002]. These instabilities may inject energy 
into the cascade at fcypi ~ 1 and uj ~ Qi. It is not possible 
yet to determine whether such energy injection is present in 
the data. The role of kinetic instabilities will be investigated 
in future studies. 

Besides the presence of a single isotropic driving scale, 
two other key assumptions were made in the construction of 
our model: 

Critical Balance: We followed Goldreich and Sridhar 
[1995] in assuming that the turbulent cascade maintains a 
state of critical balance, uj = uOni- Imposing critical balance 
and choosing a driving wavenumber kopi enables one to de- 
termine the frequency of the cascade with respect to the ion 
cyclotron frequency. Critical balance is used to approximate 
the energy flow on the (fcx,fc||) plane as a one-dimensional 
path, as shown by the solid line in Figure 1; the turbulent 
energy is, however, expected not to be restricted to this line 
but to fill the shaded region for which uj < ujni (the hydro- 
like modes in the nomenclature of Oughton et al. [2004]); 
simulations of MHD turbulence support this notion [Cho 
and Vishniac, 2000; Maron and Goldreich, 2001; Cho et al., 
2002; Cho and Lazarian, 2003; Oughton et al., 2004]. An 
extension of this cascade model to two dimensions, allow- 
ing energy to fill modes throughout the {k±,k^^) plane, is 
required to relax this assumption; we do not believe such a 
treatment would substantially alter the conclusions of this 
paper. Ultimately, the premise of critical balance will be 
evaluated by the nonlinear gyrokinetic simulations of mag- 
netized turbulence. 

Linear Damping Rates for Nonlinear Cascade: We 
assumed that the damping rates for linear eigenmodes are 
appropriate for the nonlinear turbulent cascade in which the 



energy in a wave mode is transferred to higher wavenum- 
bers within a single wave period. Li et al. [2001] call into 
question the validity of this linear damping approach, citing 
results in contrast to observation; for a more detailed com- 
parison of this paper to their work, see § 4.2. We find that 
our results are consistent with a range of observations, but 
we must leave the final judgment to nonlinear gyrokinetic 
simulations. 

Finally, we ignored a number of potentially significant 
physical effects. One of them was the possible presence 
of small-scale plasma instabilities briefiy mentioned above. 
Some of the others are discussed below. 

Radial Expansion of the Solar Wind: Not only 
does this have the effect of altering the parameters of the 
plasma — for example, the plasma transits from the low /3i 
corona to the order unity f3i solar wind at 1 AU — but the 
expansion may also affect the turbulent fluctuations. We 
focus here on the properties of the slow solar wind at 1 AU, 
but admit a possible influence of this expansion on our re- 
sults. There exists a substantial literature focusing on the 
spatial transport and evolution of the solar wind through- 
out the heliosphere [see Tu and Marsch, 1995, for a review]. 
The turbulence in these models is usually specifled locally 
by prescription. Our model is also local, so the effect of the 
expansion on the turbulence is not addressed. 

Imbalanced Turbulence: Early measurements of the 
fast solar wind at large scale suggested it was dominated by 
anti-sunward propagating Alfven waves [Belcher and Davis, 
1971]. Subsequent measurements in the fast solar wind 
found that anti-sunward Elsasser energy fluxes can exceed 
sunward fluxes by more than an order of magnitude; in 
the slow wind, the imbalance is reduced, sometimes show- 
ing equal energy in both directions [Tu et al, 1990; Tu and 
Marsch, 1990; Grappin et al., 1990; Tu and Marsch, 1995]. 
Our model only had one flux quantity e, i.e., it assumed 
a balance of Elsasser energy. An imbalanced turbulence re- 
quires a somewhat more complicated treatment [Dobrowolny 
et al, 1980; Grappin et al., 1983; Politano et al., 1989; Hos- 
sain et al., 1995; Lithwick and Goldreich, 2003; Lithwick 
et al., 2007]; thus, a more general cascade model would need 
to include an additional parameter to account for the imbal- 
ance between the fluxes. While this means that our model 
is likely not appropriate for the fast solar wind, we believe 
that it is applicable to the slow solar wind when the Elsasser 
energy fluxes are balanced. 

Note also that our model does not address the implica- 
tions of the conservation of magnetic and cross helicities in 
MHD. 

Compressive Modes: Compressive fluctuations, such 
the MHD slow, entropy, and fast modes, or, more generally, 
the density and fleld strength fluctuations, are left out of 
our model. Observational evidence suggests less than 10% of 
the power resides in compressive modes [Belcher and Davis, 
1971]. Linear theory predicts that these modes are damped 
when fcyAmfpi > 1 in the weakly coUisional solar wind [Bra- 
ginskii, 1965; Barnes, 1966]. However, it has been suggested 
that a fast wave cascade may be able to transfer energy 
to high frequencies where it can heat ions perpendicular to 
the magnetic fleld [Ghandran, 2006]. The cascade of slow 
and entropy modes is also under investigation [Schekochihin 
et al., 2007] — it is not clear that the parallel cascade of these 
modes is similar to the Alfven wave cascade. If a cascade 
of compressive fluctuations carries a significant amount of 
energy to the ion Larmor radius scale, this can contribute 
to the kinetic Alfven wave cascade in the dissipation range, 
thus altering somewhat the energetics of our model. 

Kinetic Entropy Cascade: We have referred to colli- 
sionless kinetic damping of electromagnetic fiuctuations as 
the mechanism that leads to plasma heating. The way this 
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heating occurs is, in fact, quite complex. CoUisionless damp- 
ing is an energy conversion mechanism that transfers the 
electromagnetic energy of the MHD and kinetic Alfven wave 
cascades into purely kinetic cascades of ion and electron en- 
tropy fluctuations — cascades that occur in phase space and 
that take the energy to sufficiently small spatial and velocity 
scales for the collisions to render the "coUisionless" damping 
irreversible and to effect the heating. How these cascades 
work is explained in Schekochihin et al. [2007]. An impor- 
tant point here is that these cascades should have observa- 
tional signatures, in particular, steep power-law spectra of 
low-frequency electric and magnetic fluctuations [Schekochi- 
hin et al., 2007]. A superposition of these kinetically induced 
fluctuations with the kinetic Alfven wave cascade may con- 
tribute to the wide spread in the observed spectral index in 
the dissipation range. 

4.2. Relation to Previous Work 

The cascade model proposed here is similar to that con- 
structed by Zhou and Matthaeus [1990b] and later used by 
Li et al. [2001] in a study of the dissipation of solar wind 
turbulence. The Zhou and Matthaeus [1990b] model em- 
ployed a diffusion operator for the nonlinear cascade of 
energy in wavenumber space [Leith, 1967] to develop an 
isotropic model for turbulence in the MHD regime. This 
model was used by Li et al. [2001] to determine the charac- 
teristics of the resulting energy spectrum in the dissipation 
range when coUisionless damping by the cyclotron and Lan- 
dau resonances was considered; they concluded that these 
mechanisms could only cause a steep cutoff in wavenumber 
rather than the power-law seen in observations. We come to 
a different conclusion for reasons that we shall now explain. 

First, the Li et al. [2001] model did not assume that the 
cascade frequency was determined by critical balance, so 
the change of dynamics leading to a faster cascade begin- 
ning at fcxpi ~ 1 was neglected. Second, they followed the 
cascade along constant angles in (fcx,fc||) plane of S = 0°, 
30°, 45°, and 60°, effectively assuming either isotropy or a 
scale-independent anisotropy. On the other hand, our criti- 
cally balanced model predicts that at k^pi ~ 1 — where col- 
lisionless damping first becomes significant — the anisotropy 
should be fcy/fcx ~ 0.1, leading to an angle 6 > 85°. For 
these nearly perpendicular modes, the frequency remains 
low, u! <^ Q,i, so the Alfven wave cascade is damped only 
by the Landau resonance, and the damping is less vigorous. 
Hence, we find a cutoff due to coUisionless damping that is 
much less sharp than that found by Li et al. [2001], and that 
is consistent with in situ measurements. 

In a related work, Gary and Borovsky [2004] also found 
that ion cyclotron damping produced a steep damping in- 
consistent with observations. Furthermore, they noted that, 
in the limit that fcy < fcx, electron Landau damping would 
dominate coUisionless damping with a more gradual onset. 
They suggested a qualitative picture that, in this case, the 
location of the spectral break between the inertial and dissi- 
pation ranges would be dependent on the turbulent energy 
cascade rate. We agree with their prediction of the more 
gradual onset of Landau damping, but believe the spec- 
tral break is most often caused by a change in dynamics 
at fcxPi ~ 1, a feature not included in their paper. 

Smith et al. [2001] investigated an unusual period of so- 
lar wind turbulence with jSi <^ 1, a condition in which the 
ion inertial length and the ion Larmor radius are well sep- 
arated; they concluded that the dissipation processes occur 
at the ion inertial length and that this finding is inconsis- 
tent with dissipation due to the ion cyclotron resonance. 
Here we compare our hot plasma cascade model predic- 
tions directly to interval (c) of Smith et al. [2001]. From 
Figure 1 of their paper, we estimate the following mean 



values of the parameters over interval (c): Bo — 6 nT, 
Vaw ^ 350 km/s, ~ 6 X 10* K, and fii ~ 0.01. We 
employ our fiducial driving wavenumber fcoPi = lO"**, and 
take the mean electron temperature for the slow solar wind 
as Te ~ 1.3 X 10^ K [Newbury et al, 1998], giving a tem- 
perature ratio Ti/Te = 0.5. For these parameters, the ion 
Larmor radius is pi ~ 5.5 x 10^ cm. As presented in the 
upper panel of Figure 14, the steady-state magnetic en- 
ergy spectrum from the hot plasma cascade model breaks 
from the —5/3 spectral index at a value of fcx/O; — 0.25; 
we transform from the wavenumber spectrum generated by 
our model to a frequency spectrum in the satellite frame by 
the Taylor hypothesis, using / = (fcxPi)wsm/(27rpi) to find 
a breakpoint frequency of / ~ 0.25 Hz. Smith et al. [2001] 
measure the breakpoint frequency to be / ~ 0.168 Hz; the 
hot plasma cascade model therefore predicts a breakpoint 
frequency that is roughly consistent with the observations. 
In the lower panel of Figure 14, it is shown that electron 
Landau damping is responsible for the break. In this case, 
the damping is sufficient to prevent the cascade from reach- 
ing the ion cyclotron frequency; thus, ion cyclotron damp- 
ing is negligible, falling below the plotted range of the lower 
panel. We note here that the gyrokinetic cascade model 
gives results nearly indistinguishable from those presented 
in Figure 14 for the hot plasma cascade model; this is not 
surprising since the cyclotron resonance plays a negligible 
role. It is worth noting that the wavenumber associated 
with the breakpoint, when cast in terms of the ion inertial 
length, is fcxdi — 2.5, so we agree with the Smith et al. 
[2001] conclusion that the breakpoint occurs at the ion in- 
ertial length scale. The salient point here is that, although 
the magnetic energy spectrum breakpoint often occurs at 
fcxPi ~ 1 due to the change in linear physics from the MHD 
to the kinetic Alfven wave regime, strong damping may shift 
the breakpoint to a lower value of k±pi. 

4.3. Application to the Solar Corona 

Although this paper has focused on anisotropic turbu- 
lence in the solar wind, such turbulence is also expected 
to be present in the solar corona. One of the key prob- 
lems for the solar corona is the identification of the mecha- 
nism by which the coronal plasma is heated to yield the ob- 
served minority ion temperatures higher than proton tem- 
peratures T|qji > Tp and large ion and proton tempera- 
ture anisotropics Tx > Ty [Kohl et al, 1997, 1998; An- 
tonucci et al, 2000; Kohl et al, 2006]. In the spirit of a 
long line of previous investigations [Hollweg, 1986; Isenberg, 
1990; Li et al, 1999; Matthaeus et al, 1999; Dmitruk et al, 
2001, 2002; Cranmer and van Ballegooijen, 2003, 2005; 
Cranmer et al, 2007], we brifely address what our cascade 
model predicts for ion heating in the solar corona. 

The first caveat is that the cascade model described in 
this paper assumes a balance of the oppositely-directed 
Elsasser energy fiuxes along the magnetic field line. There- 
fore, the predictions of this simple model cannot be applied 
to the coronal hole regions in which the field lines are open 
and the anti-sunward Alfven wave energy fiux dominates; 
the conclusions here will only apply to regions of closed field 
lines in which the turbulence demonstrates balanced energy 
fiuxes. For the solar corona at a heliocentric radius r = 2Rq , 
we take the following plasma parameters: magnetic field 
magnitude Bq = 0.8 G, number density n = 3.7 x 10® cm~^, 
electron temperature Te = 8 x 10® K, and proton tempera- 
ture Te = 2 X 10® K. These plasma conditions give an ion 
Larmor radius pi = 2.4 x 10"^ cm, an ion cyclotron frequency 
Q.i = 7.7 X 10^ rad/s, an ion plasma beta (Hi = 0.004, and a 
temperature ratio Ti/T^ = 2.5. 

In addition to these parameters, for the hot plasma cas- 
cade model we must also specify the driving by supplying 
the characteristic isotropic driving wavenumber fcoPi. From 
the model of Cranmer and van Ballegooijen [2005], we find 
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that, at a radius of r — 2Rq, the driving energy due to pho- 
tospheric motions peaks at a period of P = 120 s. Using the 
dispersion relation for Alfven waves uj = fcyVA, this gives a 
parallel driving length scale I/y = 3.4 x 10^° cm. At this 
radius, the Cranmer and van Ballegooijen [2005] model es- 
timates the perpendicular driving scale L±_ = 1.3 x 10® cm 
(using their best value of A = 0.35 in their equation [51]); 
for the period P — 120 s, this implies a perpendicular ve- 
locity of — 110 km/s. The effective isotropic driving 
wavenumber can be found using equation (13), giving the 
result kopi = 7 x 10^^". Note that, although the length 
scale associated with fco is Lq = 2 x 10^^ cm, there are no 
motions at this scale; the driving occurs anisotropically at 
L|| > L± and the value of koPi calculated here is simply a 
convenient way to specify the critical balance relation deter- 
mined by the driving mechanism, as explained in § 1.1. 

The heating resulting from the damping of the turbulence 
can occur via either the Landau or cyclotron resonances in 
the hot plasma cascade model; Landau damping leads to 
an increase of the parallel temperature, cyclotron damp- 
ing to an increase of the perpendicular temperature. For 
the low Pi conditions in the solar corona. Landau damping 
onto the ions is negligible, so all heating is due to either 
electron Landau damping or ion cyclotron damping. Using 
equation (15) to estimate the frequency of the cascade at 
k±pi — 1 (and neglecting damping by taking e — eo), we 
find uj/fli = 1.7 X 10"^; this suggests ion cyclotron damping 
will be negligible in this case. The steady-state solution of 
the hot plasma cascade model for these parameters confirms 
that the Alfven wave cascade driven by photospheric mo- 
tions results in only parallel electron heating via the Landau 
resonance. In fact, for all values of kopi < 10~^, the turbu- 
lent cascade fails to reach the ion cyclotron frequency and 
therefore no ion heating occurs. Significant perpendicular 
ion heating via the balanced, Alfven wave cascade requires 
an isotropic driving scale Lo < 1.5 x 10^ cm, a condition 
not satisfied in the solar corona. In conclusion, this turbu- 
lent cascade model cannot account for the perpendicular ion 
heating observed in the solar corona. 

5. Conclusions 

In the study of solar wind turbulence, a central issue is to 
understand the mechanisms by which the energy contained 
in large-scale turbulent motions is transferred to smaller 
scales and higher frequencies where it is eventually deposited 
as thermal energy in the plasma. In § 2.1 we develop a sim- 
ple model that traces the anisotropic cascade of Alfven wave 
energy from the MHD to the kinetic Alfven wave regime. 
The model is based on the assumption of locality of energy 
transfer in wavenumber space, the conjecture of a critical 
balance between the nonlinear interaction and linear prop- 
agation timescales, and the premise that the linear kinetic 
damping rates determine heating rates even in the presence 
of the nonlinear cascade. It provides a guide for the inter- 
pretation of observational data in the balanced, slow solar 
wind and will be useful in connecting these observations to 
nonlinear gyrokinetic simulations of turbulence. 

The anisotropic cascade model makes several predictions 
about the turbulent magnetic energy spectrum measured in 
the balanced, slow solar wind — these are presented in § 2.2- 
§ 2.5. The cascade is weakly damped in the Alfven wave 
regime k±pi <C 1 and the spectrum has the familiar k'^^"^ 
slope. In the dissipation range, at wavenumbers above the 
observed break, the magnetic energy spectrum undergoes a 
slow exponential cut-off. We show that the measured power- 
law behavior of the spectrum in the dissipation range may 
be an artifact of limited magnetometer sensitivity (except 
in the case of very weak damping) and that over the typi- 
cal range of parameters Pi and Ti/T^ in the solar wind, the 
varying strength of Landau damping naturally reproduces 



the observed variation of the effective spectral index in the 
dissipation range from —7/3 to —4 (§ 2.3-§ 2.4). The spec- 
tral break between the inertial range and dissipation range 
power laws occurs at the ion Larmor radius k±pi ~ 1 in 
plasmas with weak damping due to the dispersive nature of 
the kinetic Alfven wave; in the presence of strong damping, 
it may occur at k^pi < 1 due to damping by the Landau 
resonance and/or the ion cyclotron resonance. See Figure 2 
for a plot of the damping rate as a function of plasma pa- 
rameters at fcxPi = 1. 

In § 3 we examine the applicability of gyrokinetics to tur- 
bulence in the balanced, slow solar wind. At scales some- 
what below the driving scale, the anisotropy of the turbu- 
lence, fc|| < fcx, is supported by observations [Dasso et ai, 
2005] in the slow solar wind. Given the observed anisotropy, 
we find that a gyrokinetic description becomes valid from 
a scale a couple of orders of magnitude below the driving 
scale up to the point where the fiuctuation frequency ap- 
proaches the cyclotron frequency, ~ Sli . The driving scale 
of the turbulence in the slow solar wind is estimated to be 
L ~ 2 X 10^^ cm (§ 3.1), in agreement with three other 
independent determinations. For typical solar wind param- 
eters, we find that fcy <^ k± and u> <^ Qi a,t a, perpendicu- 
lar wavenumber k±pi ~ 1 (§ 3.2); therefore, this transition 
regime is well described by gyrokinetics. Wave phase ve- 
locity evidence from in situ observations [Bale et al., 2005] 
argues for the existence of a kinetic Alfven wave cascade 
in the dissipation range (which is captured in our model); 
the measurements are inconsistent with the onset of ion cy- 
clotron damping (§ 3.4). 

The gyrokinetic cascade and the resultant particle heat- 
ing through Landau damping, however, cannot be the whole 
story. The selective perpendicular heating of minor ions and 
the breaking of adiabatic invariance for ions implied by ob- 
servational data, particularly in the fast wind, indicate that 
fluctuations at or above the cyclotron frequency are present 
in the solar corona and solar wind. Although a mechanism 
not incorporated into our model may be responsible — for ex- 
ample, heating by fast waves [Chandran, 2006] or velocity- 
space instabilities that transfer energy from the parallel to 
the perpendicular direction [Scarf et ai, 1967; Eviatar and 
Schulz, 1970; Kasper et al, 2002; Marsch, 2006] — we may 
employ our model to determine under what conditions the 
gyrokinetic cascade reaches the ion cyclotron frequency. For 
plasma parameters typical of the solar wind and the driv- 
ing scale estimate L ~ 2 x 10^^ cm, the gyrokinetic cas- 
cade model generally reaches the ion cyclotron frequency at 
k±pi ^ 1, by which point Landau damping has removed 
much of the cascade energy. In addition, damping of the 
cascade can slow the increase of the frequency with perpen- 
dicular wavenumber, possibly preventing the cyclotron fre- 
quency from being reached at all, as shown in Figure 8. A 
strong constraint on the wavenumber threshold at which the 
cyclotron resonance becomes important — the same thresh- 
old at which the gyrokinetic approximation is violated — 
is determined by a comparison of the gyrokinetic and hot 
plasma cascade models. The regime of validity for gyroki- 
netics is given by the shaded region in Figure 9 and cov- 
ers most of the parameter range of interest in the study 
of turbulence in the solar wind. In the range of parame- 
ters relevant to the solar wind, this threshold roughly scales 
as k±pi ~ {koPi)~^^'''Pi^^ (§ 3.5). Therefore, the cyclotron 
resonance can have a significant effect on the dissipation- 
range dynamics only in plasmas with Pi < 0.01 or in sys- 
tems driven isotropically at scales L <^ 2 x 10^^ cm. If 
the cyclotron resonance is reached within the dissipation 
range before the electron Landau damping becomes domi- 
nant, there is a conversion between kinetic Alfven waves and 
ion Bernstein waves in the narrow frequency band around 
the resonance. However, the amount of cyclotron damping 
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associated with this conversion is quite small and the nar- 
rowness of the resonance suggests that the kinetic Alfven 
wave cascade may in fact continue virtually unimpeded even 
at a; > rJi (see § 3.6 and Appendix A). 

In § 4.3, we apply this cascade model to the problem of 
ion heating in the solar corona. For the low /3i conditions 
of the corona, the dissipation of turbulence can cause par- 
allel heating of the electrons via the Landau resonance or 
perpendicular heating of the ions via the ion cyclotron res- 
onance. We conclude, however, that the highly anisotropic 
driving due to photospheric motions results in only parallel 
electron heating via the Landau resonance; only for an ef- 
fective isotropic driving scale Lo < 1.5 x 10^ cm, a condition 
not satisfied in the corona, will the Alfven wave cascade as 
modeled here lead to substantial perpendicular ion heating. 

Although our simple cascade model allows us to produce 
plausible predictions for magnetic-energy spectra, reaching a 
definitive understanding of turbulence in the solar wind will 
require nonlinear simulations of magnetized plasma turbu- 
lence in the weakly colhsional regime. A program of gyroki- 
netic simulations is currently underway [Howes et al, 2008; 
Howes, 2008]. These simulations will both test the validity 
of the ideas that underpin the approach taken in this paper 
and probe the physics beyond the scope of our simple model, 
in particular the complicated dynamics at the scale of the ion 
Larmor radius fcxPi ~ 1- To provide detailed information 
for comparison to in situ measurements, gyrokinetic sim- 
ulations will calculate particle heating rates, magnetic and 
cross helicities, and nonlinear energy transfer rates in imbal- 
anced cascades. The cascade model proposed in this paper 
is a first step in connecting the results of nonlinear gyroki- 
netic simulations to the observations in the solar wind and 
in assessing the applicability of gyrokinetics to this problem. 
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Appendix A: Plasma Dispersion Relation, 
Mode Conversion and Damping at the Ion 
Cyclotron Resonance 

The general eigenvalue problem for linear perturbations 
of a thermal plasma is 



D ■ (5E = 0. 



(Al) 



where, if we choose the coordinate frame so that z is the 
direction of the mean magnetic field and x is the direction 
of kx, the components of the dielectric tensor are (see, e.g., 
Stix 1992) 
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where Qs — k^Ps, Cns ~ {lu — nQs)/\k\\\vths, In are modified 
Bessel functions, and Z is the plasma dispersion function. 

We now consider the general hot plasma dispersion 
detD = 0, or 

DxxDyy DzZ H~ '^DxyDyzDxZ 

- Dl.Dyy + Dl.D^x + DlyD,, = 0, (A8) 

in the asymptotic limits that we believe to be appropri- 
ate for the dissipation range of the solar wind: k±pi ^ 1, 
k±pe <^ 1, and fcy <^ k±. Then we may expand the Bessel 
and plasma dispersion functions in the above expressions for 
Dij, taking fcy < k±, Oe < 1, Coe < 1, Cie » 1 for n / 0, 
Qi 3> 1, Cm 3> 1. After straightforward algebra, we get 



= [1 - 2^' (1 - C)] , 
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where ui = u!/\k^^\vAk_Lpi. The (small) terms containing 

(A15) 
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are responsible for electron Landau damping of the kinetic 
Alfven waves, while the effect of the cyclotron resonance is 
(to lowest order) contained only in the matrix element D^^ 
(corresponding to perpendicular electrostatic fluctuations): 
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where Cli = Qi / {\k^^VA\k±pi) and, in terms of the normal- 
ized frequencies uj and Qi, C,ni = (Co — nQ.)k±pi/ ^fpl. The 
dispersion relation (A8) can then be written as follows: 
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If we set C — here, we recover the gyrokinetic dispersion 
relation (in the limit k±pi ';$> 1, k±pe <^ 1), whose solutions 
are the weakly Landau-damped kinetic Alfven waves (see 
section 2.6.2 of Howes et al. [2006]). The normalized fre- 
quency and damping rate of these waves in our notation are 
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The kinetic Alfven waves will be affected by the ion cy- 
clotron resonance when their frequency is close the cyclotron 
frequency, i.e., Qi ~ ljkaw, or 



-{k±pi) 



(A20) 



For the highly anisotropic fluctuations in the kinetic Alfven 
wave cascade, we have fcy/fcx <^ 1, so condition (A20) is 
reached deep in the dissipation range, k±pi ^ 1 — where ex- 
actly depends on the degree of anisotropy of the fluctuations 
and can be estimated using the critical balance assumption 
to relate fcy and k± (see equation (16), which also includes 
the dependence on ft and Ti/Te). Note that gyrokinetics 
breaks down at this point. 

Let us consider what happens in the region of the 
wavenumber space where the condition (A20) is satisfied. 
Define A„ = wkaw — nCli and assume A„ ^ nCl ~ iDkaw ~ 
1. We now look for solutions of the dispersion relation (A17) 
in the form Co — cjkaw + where 5Cj ~ A„ ^ i^kaw (the 
validity of this assumption will be confirmed shortly) . In the 
dispersion relation (Al7), we neglect the Landau damping 
terms (they are controlled by an independent small param- 
eter, L oc \J nie/mi) and get 
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We shall see that the imaginary part of C is exponentially 
small, so to lowest order, equations (A21) and (A22) com- 
bine into a quadratic equation for Suj. The solutions of this 
equation are 
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When the kinetic Alfven wave frequency is not too close 
to the cyclotron resonance, namely |An| 3> 1/^/kTpi, the 
"+" solution corresponds to the kinetic Alfven wave 
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and the "— " solution to an approximately electrostatic wave 
called the ion Bernstein wave (see, e.g., Stix 1992), whose 
frequency is close to nQ,i: 
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(this solution is high frequency and is not captured by gy- 
rokinetics) . 

When jA„| ~ l/^k±pi, the more general formula (A23) 
describes the conversion of the kinetic Alfven wave into the 
ion Bernstein wave below the cyclotron resonance (iDkaw < 



fii) and of the ion Bernstein wave into the kinetic Alfven 
wave above the resonance (wkaw > f^O- This transition 
is illustrated in Figure 10 for the n = 1 resonance. The 
behavior at n > 1 is similar. 

The damping rate can now be calculated perturbatively; 
restoring the Landau term L and the imaginary part of C 
and replacing Sub ^ Sub -\- 47 in equations (A21) and (A22) 
with 5u) given by equation (A23), we get (keeping only the 
lowest order terms for k^pi ^ 1) 
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where 7kaw is the electron Landau damping rate given 
by equation (A 19) an d ("„i — (A„ + 5Cj)k ±pi/ ^ /JSj. Since 
Slu ~ A„i ~ l/v^fcj~pi, we find that C^ni ~ x/fcjTpT, so the cy- 
clotron damping remains exponentially small and we have 
thus confirmed a posteriori that the expansion of Z{C^ni) in 
large argument was justified. The electron Landau damping 
rate (the first term in equation (A26)) and the ion cyclotron 
damping rate (the second term in equation (A26)) are plot- 
ted in Figure 10. Note that the ion cyclotron damping is 
quite weak both in absolute terms and in comparison with 
the electron Landau damping. 

The main points that we have demonstrated by the above 
calculation are itemized at the beginning of § 3.6, where they 
are followed by a discussion of a somewhat less asymptotic 
case. That discussion and a much more detailed numeri- 
cal study by Li and Habbal [2001] demonstrate some of the 
possibilities for mode conversion and damping at the ion 
cyclotron resonances that are beyond the scope of a simple 
asymptotic calculation like the one presented in this Ap- 
pendix. 

Notes 

1. Note that, although we parameterize the cascade in terms of 
the isotropic driving wavenumber fcoPi, this parameterization 
does not require that the cascade be driven isotropically. A cas- 
cade driven both in critical balance and anisotropically, with 
^-L 7^ ''^ll I can be characterized by an effective value of fcopi 
along the extrapolated cascade path at which fcx = k.\i^. 

2. The Vlasov-Maxwell dispersion relation for a hot plasma agrees 
with the gyrokinetic results in the limit fcy <C fcx and w ^ O;, 
precisely the anisotropic and low-frequency fluctuations antic- 
ipated by the GS theory for MHD turbulence. 

3. The Landau resonance includes both Landau damping by the 
parallel electric field perturbation and transit-time, or Barnes 
[1966] , damping by the parallel magnetic field perturbation; in 
this paper, we refer to the collective damping effect of these 
mechanisms as Landau damping. 

4. There is a typo in the legend of this figure: the spectrum at 
0.9 AU is for the slow solar wind (L), not the fast (H). 
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Figure 1. Schematic diagram of the nonlinear transfer 
of energy across the {k±, fey) plane from the low driving 
wavenumber to high wavenumbers where the turbulence 
is dissipated by kinetic processes. Driving is chosen to be 
isotropic at a wavenumber kopi — 10~^. Critical balance, 
ui = Urii, defines a one-dimensional path for the turbulent 
cascade of energy on the (fcx,fc||) plane denoted by the 
solid line. Numerical simulations show turbulent energy 
over the entire shaded region, in which ui < uini- The 
cascade model constructed in this paper can be thought 
of as describing energy integrated over all fc|| at each value 
of k±. 
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Figure 2. Contours of log(7/ti;) on the plane {pi, Ti/Te) 
for k^pi = 1. The solid shading corresponds to the 
range of parameters in the slow solar wind, the dashed 
shading to the fast solar wind. Contours are spaced at 
Alog(7/ci;) = 0.1 with thick dashed lines at each unit of 
the logarithm; contour values for each logarithmic unit 
are denoted. 
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Figure 3. One-dimensional magnetic energy spectra 
for tlie turbulent cascade from the MHD to the kinetic 
Alfven wave regime. The chosen parameters are j3i = 1 
and Ti/Te — 1. Solutions of our cascade model with no 
damping (dashed) and with the gyrokinetic damping rate 
(solid) are shown. 
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Figure 4. One-dimensional magnetic energy spectra for 
three gyrokinetic models: (1) Pi = 0.5, Ti/T^ = 3, (2) 
Pi = 3, Ti/n = 0.6, (3) ft = 0.03, T./Te = 0.175. Panel 
(a) shows that all three spectra demonstrate a dissipa- 
tive roll-off that cannot be fit by a single spectral index 
k±pi > 1 (consistent with an exponential roll-off). 
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Figure 5. One-dimensional magnetic energy spectra for 
three gyrokinetic models: (1) /3, = 0.5, Ti/Te = 3, (2) 
Pi = 3, Ti/n = 0.6, (3) ft = 0.03, T./Te = 0.175. 
Panel (b) shows the same spectra as panel (a) with an 
added constant sensitivity limit, which causes resulting 
the spectra to resemble power laws. The sensitivity limit 
is two orders of magnitude below the value of the spec- 
trum at k±pi = 1 for spectrum 1, three for spectra 2 
and 3. 
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Figure 6. The effect of changing the Kolmogorov con- 
stants over a range of values C\^^C2 ~ 3, 6, and 12. The 
measured effective spectral index in the dissipation range 
varies by +1/3 or— 2/3. 
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Figure 7. The normalized wave phase velocity uj/{k\\VA) 
vs. wavenumber kpi for models with different assump- 
tions about the anisotropy of the fluctuations. The cal- 
culations use the full linear plasma dispersion relation 
for a coUisionless plasma with Pi = 1 and Ti/T^ — 
1. We assume a critically balanced cascade with the 
parallel wavenumber governed by equation (11) with 
no dissipation, e = eo, and driving scales kopi = 
10~^, 10^'', lO^'^, and 1. The rise in wave phase veloc- 
ity at kpi ~ 1 for kopi <^ 1 is due to the transition to 
kinetic Alfven waves; this is consistent with the in situ 
measurements by [Bale et al, 2005]. The turn-down in 
wave phase velocity, due to the onset of ion cyclotron 
damping, is not observed. 
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Figure 8. The frequency tj/Sli [equation (15)] and par- 
allel wavenumber fcypi vs. k±pi [equation (11)] for un- 
damped (dashed) and gyrokinetic (solid) steady-state so- 
lutions of the gyrokinetic cascade model. If the damping 
of nonlinear turbulent fluctuations is less than or equal 
to the linear gyrokinetic damping rates, than the non- 
linear frequency should lie within the vertically shaded 
region and the parallel wavenumber within the diago- 
nally shaded region. The effective damping clearly plays 
a strong role in determining whether the turbulent fluc- 
tuations reach the cyclotron frequency, ui/fli ~ 1. 
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Figure 9. The threshold value of k_i_pi at which the 
ion cyclotron resonance has an effect on the linear wave 
modes with fcopi = 10^'' and Ti/Te = 1- The threshold 
is determined by comparing the solutions of the gyroki- 
netic and hot plasma dispersion relations. The solid and 
dashed lines correspond to the path in the {k±, fc||) plane 
calculated from the gyrokinetic model with and without 
damping, respectively; the dotted line is the analytical es- 
timate given by equation (16) (only plotted for the typical 
range of j3i in the solar wind near Earth). 
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Figure 10. Approximate solution of the hot plasma dis- 
persion relation illustrating the conversion between ki- 
netic Alfven and ion Bernstein waves (see Appendix A). 
Plotted in the upper panel are frequencies uj/Qi [equa- 
tion (A23)] of the low- (solid) and high-frequency 
(dashed) branches. The characteristic wave modes along 
each branch are labeled as kinetic Alfven wave (KAW) 
and ion Bernstein wave (IBW). Shown in the lower panel 
are the electron Landau damping (electron) 7e/r2i and 
ion cyclotron damping (ion) 7^ /f2i , given by the first and 
second terms in equation (A26), respectively. The plasma 
parameters are /3i — 0.3 and Ti/Te = 1. The solutions are 
plotted along the critically balanced cascade path with 
k\\{k±) given by equation (11), where fcopi = 10"'' and 
uj = a is given by the second expression in equation (4). 
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Figure 11. Real frequencies from the solutions of the 
hnear gyrokinetic dispersion relation (dashed line) and 
the hot plasma dispersion relation (dotted line) for /3i — 
0.03, Ti/Te = 3, and fcopi = 2 x 10"*. The solid line 
shows the path in the (fcx,fc||) plane implied by critical 
balance. Two branches of the hot plasma solution are 
shown, corresponding, in the k^pi limit, to the MHD 
Alfven wave and the MHD fast wave. 
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Figure 12. Normalized ion and electron damping rates, 
7i/a; and "fs/i^, from the solutions of the linear gyroki- 
netic dispersion relation (dashed line) and the hot plasma 
dispersion relation (dotted line). In the upper panel we 
show the low-frequency (Alfven) branch solutions; in the 
lower panel, high-frequency (fast) branch. The gyroki- 
netic solution shows only the electron damping — the ion 
damping rates are below the plotted range. For the hot 
plasma solutions, the contributions from ion and electron 
damping are labeled accordingly. These solutions are for 
Pi = 0.03, Ti/Te = 3, and kop, = 2 x 10"" and the 
cascade path in the {k±, fcy) plane shown in Figure 11. 
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Figure 13. One-dimensional magnetic energy spectra 
for /3i = 0.03, Ti/Te = 3, and kop^ = 2 x 10""*. Gy- 
rokinetic and hot plasma cascade models are compared; 
both models show similar spectra in the kinetic regime 
with the effective spectral index in the dissipation range 
of approximately —4. 
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Figure 14. The upper panel presents the normalized 
one-dimensional magnetic energy spectrum using the hot 
plasma cascade model for /3i = 0.01, Ti/Te = 0.5, and 
kopi = 10~*, the parameters suitable for modeling the 
interval of solar wind studied by Smith et al, [2001]. The 

"5/3 

spectrum breaks from the behavior at fcxPi — 0.25, 

as indicated by the arrow. The lower panel shows the fre- 
quency of the turbulent fluctuations w/Oi and the nor- 
malized electron damping rate 7e /a; vs. k±pi. The damp- 
ing is sufficient to prevent the cascade from reaching the 
ion cyclotron frequency SI;; therefore, the break in the 
spectrum at k±pi ~ 0.25 is due to damping via the Lan- 
dau resonance with the electrons. 
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Table 1. Definitions 

Species s 

Speed of light c 

Mass TTla 

Charge Qs 

Number Density Us 

Temperature'' Ts 
Species Plasma Beta ~ StvUsTs / Bq 



Plasma Frequency cups = ^TiUsqi / ms 

Cyclotron Frequency Sis = qsBo/{msc) 

Thermal Velocity Vtha = \/2Ts/ms 

Alfven Velocity va = Bo^/Aivnimi 

Mean Free Path A^fpa 

Larmor Radius ps = vths/^s 

Inertial Length ds — c/ujps 

Temperature is given in energy units, absorbing the Boltzmann constant. 



